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Abstract 

Recent experimental and theoretical advances in investigating electromag- 
netic meson production reactions in the nucleon resonance region are reviewed. 
The article gives a description of current experimental facilities with electron 
and photon beams and presents an unified derivation of most of the phe- 
nomenological approaches being used to extract the resonance parameters 
from the data. The analyses of ir and rj production data and the result- 
ing transition form factors for the A(1232)P 3 3, N(1535)Sii, N(1440)Pn, and 
N(1520)Di3 resonances are discussed in detail. The status of our understand- 
ing of the reactions with production of two pions, kaons, and vector mesons 
is also reviewed. 



1. Introduction 

The quest for understanding the structure and interaction of hadrons has been 
the motivation of strong interaction physics for decades. The advent of Quan- 
tum Chromodynamics (QCD) 1 has led to a general theoretical description of the 
strong interaction in terms of the fundamental constituents, quarks and gluons. At 
very high energies, perturbative methods have proven very effective in the descrip- 
tion of many processes. However, because of the complexity of the theory, we are 
still a long way from being able to describe the strong force as it is manifest in 
the structure of baryons and mesons. The most fundamental approach to resolve 
this difficulty is to develop accurate numerical simulations of QCD on the Lattice 
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(Lattice QCD) 2 . Alternatively, hadron models with effective degrees of freedom 
have been constructed for interpreting data. For example, near threshold pion-pion 
scattering, pion-nucleon scattering, and pion photoproduction can be successfully 
described by chiral perturbation theory 3 which is formulated in terms of hadron 
degrees of freedom and constrained only by the symmetry properties of QCD. The 
constituent quark model 4,5 is another successful, though not fully understood, ex- 
ample. In some cases, the results from these two different theoretical efforts arc 
complementary in understanding the data and making predictions for future exper- 
iments. 

For heavy quark systems, Lattice QCD (LQCD) can now predict accurate quan- 
tities for interpreting the data from, for example, B meson facilities. For light-quark 
systems, the small quark masses are difficult to implement, and approximations have 
been necessary in Lattice QCD calculations. Nevertheless, significant progress has 
been made in calculating some basic properties of baryons, such as masses of ground 
states, as well as of low lying excited states 6,7,8,9 . Even the first LQCD calcula- 
tion of the electromagnetic transition form factors from the ground state proton to 
the first excited state, the A(1232), has been attempted recently 10 . However, reli- 
able Lattice QCD calculations for electromagnetic meson production reactions, the 
subject of this article, seem to be in the distant future. In the foreseeable future, 
models of hadron structure and reactions will likely continue to play an important 
role and provide theoretical guidance for experimenters. 

The development of hadron models for the nuclcon and nucleon resonances 
(N*) has a long history. In the past three decades, the constituent quark model 
has been greatly refined to account for residual quark-quark interactions due to 
one-gluon-exchange 4 ' 5 and/or Goldstone boson exchange 11,12 . Efforts are under- 
way to re- formulate the model within the relativistic quantum mechanics 13,14,15 . 
Conceptually completely different models have also been developed, such as bag 
models 16 , chiral bag models 17,18 , algebraic models 19 , soliton models 20 , color dielec- 
tric models 21 Skyrme models 22 , and covariant models based on Dyson-Schwinger 
equations 23 . With suitable phenomenlogical procedures, most of these models are 
comparable in reproducing the low-lying N* spectra as determined by the amplitude 
analyses of elastic irN scattering. However they have rather different predictions 
on the number and ordering of the highly excited N* states. They also differ sig- 
nificantly in predicting some dynamical quantities such as the electromagnetic and 
mesonic N-N* transition form factors. Clearly, accurate experimental information 
for these N* observables is needed to distinguish these models. This information 
can be extracted from the data of electromagnetic meson production reactions. In 
the past few years, such data with high precision have been extensively accumulated 
at Thomas Jefferson National Accelerator Facility (JLab), MIT-Bates, and LEGS of 
Brookhaven National Laboratory in the United States, MAMI of Mainz and ELSA 
of Bonn in Germany, GRAAL of Grenoble in France, and LEPS of Spring-8 in 
Japan. In this paper we will review these xperimental developments and the status 
of our understanding of the data accumulated in recent years. Our focus will be on 
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the study of N* excitations. The use of these data for other investigations will not 
be covered. 

It is useful to briefly describe here the recent advances in using the new data to 
address some of the long-standing problems in the study of N* physics. The first 
one is the so-called missing resonance problem. This problem originated from the 
observation that some of the N* states predicted by the constituent quark model 
are not seen in the baryon spectra determined mainly from the amplitude analyses 
of nN elastic scattering. There are two possible solutions for this problem. First, it 
is possible that the constituent quark model has wrong effective degrees of freedom 
of QCD in describing the highly excited baryon states. Other models with fewer 
degrees of freedom, such as quark-diquark models or models based on alternative 
symmetry schemes 206 , could be more accurate in reproducing the baryon spectra. 
The second possibility is that these missing resonances do not couple strongly with 
the ttN channel and can only be observed in other processes, as suggested by Isgur 
and Koniuk 169 in 1980. Data from the experiments measuring as many meson- 
baryon channels as possible are needed to resolve the missing resonance problem. 

The second outstanding problem in the study of N* physics is that the partial 
decay widths of baryon resonances compiled and published periodically by the Par- 
ticle Data Group (PDG) have very large uncertainties in most cases 24 . For some 
decay channels, such as rjN, KT, and ojN , the large uncertainties are mainly due 
to insufficient data. But the discrepancies between the results from using differ- 
ent amplitude analysis methods is also a source of the uncertainties. This problem 
can be resolved only with a sufficiently large data base that allows much stronger 
constraints on amplitude analyses, and a strong reduction of the model dependence 
of the extracted partial decay widths as well as other N* parameters. This re- 
quires that the data must be precise and must cover very large kinematic regions 
in scattering angles, energies, and momentum transfers. The data of polarization 
observables must also be as extensive as possible. 

The above two experimental challenges have been met with the operations of the 
electron and photon facilities mentioned above. These facilities are also equipped 
with sophisticated detectors for measuring not only the dominant single pion chan- 
nel but also kaon, vector meson, and two-pion channels. The CEBAF Large Accep- 
tance Spectrometer (CLAS) at JLab is the most complete and advanced detector 
in the field. 

The third long-standing problem is in the theoretical interpretations of the N* 
parameters listed by the PDG. Most of the model predictions on N* — > 77V helicity 
are only in a very qualitative agreement with the PDG values. In some cases, they 
disagree even in signs. One could attribute this to the large experimental uncer- 
tainties, as discussed above. However, the well determined empirical values of the 
simplest and most unambiguous A — > helicity amplitudes are about 40% larger 
than the predictions from practically all of the hadron models mentioned above. 
This raises the question about how the hadron models as well as the Lattice QCD 
calculations are related to the N* parameters extracted from empirical amplitude 
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analyses. We need to evaluate critically their relationships from the point of view 
of fundamental reaction theory. The discrepancies in the A(1232) region must be 
understood before meaningful comparisons between theoretical predictions and em- 
pirical values can be made. Much progress has been made in this area. The results, 
as will be detailed in section 5.1, strongly indicate that it is necessary to apply 
an appropriate reaction theory in making meaningful comparisons of the empirical 
values from amplitude analyses and the predictions from hadron models and LQCD. 

Summing up the above discussions, it is clear that in the absence of a funda- 
mental solution of QCD in the resonance regions, the study of N* excitations needs 
close collaborations between theoretical and experimental efforts. This is illustrated 
in Fig. 1. On the theoretical side, we need to use Lattice QCD calculations and/or 
hadron structure models to predict properties of nucleon resonances, such as the 
N-N* transition form factors indicated in Fig. 1. On the experimental side, we 
need to accumulate sufficiently extensive and precise data of meson production re- 
actions. We then must develop reaction models for interpreting the data in terms 
of hadron structure calculations. The development of empirical amplitude analyses 
of the data is an important part of this task. 



In section 2, the current experimental facilities will be reviewed. The general 
formulation for calculating cross sections of electromagnetic meson production is 
presented in section 3. Section 4 is devoted to provide a unified derivation of 
models used in the interpretation of the data. Results are presented in section 5. 
Concluding remarks and outlook are given in section 6. 




Figure 1: Scheme for N* study 
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2. Experimental Facilities 

2.1. Thomas Jefferson National Accelerator Facility 

The Thomas Jefferson National Accelerator Facility (JLab) in Newport News, Vir- 
ginia, operates a CW electron accelerator with energies in the range up to 6 GeV 27 . 
Three experimental Halls receive highly polarized electron beams with the same, or 
with different but correlated energies, simultaneously. Beam currents in the range 
from 0.1 nA to 150 ^A can be delivered to the experiments, simultaneously. 

2.1.1. Experimental Hall A - HRS* 2 

Hall A houses a pair of identical focussing magnetic spectrometers 28 with high 
resolution (HRS 2 ), each with a momentum resolution of Ap/p ~ 2 x 10~ 4 , one 
of them is instrumented with a gas Cerenkov counter and a shower counter for 
the identification of electrons. The hadron arm is instrumented with a proton 
recoil polarimctcr. The detector package allows identification of charged pions, 
kaons, and protons. A polarized 3 He target is used for experiments that require 
polarized neutron targets. The HRS 2 spectrometers have been used to measure the 
reaction ep — > epir° in the A(1232) region and to extract various single and double 
polarization response functions. 

2.1.2. Experimental Hall B - CLAS 

Hall B houses the CEBAF Large Acceptance Spectrometer (CLAS) detector, 
and a photon energy tagging facility 29 . CLAS can be operated with electron beams 
and with energy tagged photon beams. The photon beam can be either unpolar- 
ized or can be linearly or circularly polarized. The detector system was designed 
specifically with the detection of multiple particle final states in mind. The driving 
motivation for the construction of CLAS was the nucleon resonances (N*) program, 
with the emphasis on the study of the jNN* and jNA* transition form factors, 
and the search for missing resonances. Figure 2 shows the CLAS detector. At the 
core of the detector is a toroidal magnet consisting of six superconducting coils sym- 
metrically arranged around the beam line. Each of the six sectors is instrumented 
as an independent spectrometer with 34 layers of tracking chambers allowing for the 
full reconstruction of the charged particle 3-momentum vectors. Charged hadron 
identification is accomplished by combining momentum and time-of- flight, and the 
measured path length from the target to the plastic scintillation counters which 
surround the entire tracking chambers. Timing resolutions of AT = 120 — 200 pscc 
(rms) are achieved, depending on the length of the scintillator bar which ranges 
from 30 cm to 350 cm. Mass and charge number (Z) reconstruction is shown in the 
left panel of Fig. 3. Protons and pions can be separated for momenta up to 4 GeV/c, 
and pions and kaons up to about 2 GeV/c. The wide range of particle identification 
allows to study the complete range of reactions relevant to the N* program. In 
the polar angle range of up to 70° photons and neutrons can be detected using the 
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electromagnetic calorimeters. The forward angular range from about 10° to 50° is 
instrumented with gas Cerenkov counters for the identification of electrons. 



Large-angle Calorimeter 



Drift Chambers I Electromagnetie Calorimeter 




Figure 2: The CLAS detector at JLab. Left panel: The longitudinal cut along the beam 
line shows the 3 drift chamber regions, the Cherenkov counter at forward angles, the time- 
of-flight (TOF) system, and the electromagnetic calorimeters. A simulated event shows 
an electron (upper) and a positively charged hadron. Right panel: Transverse cut through 
CLAS. The six superconducting coils provide a six sector structure with independent 
detectors. A polarized target can be inserted into the large bore near the center. 

In the N* program, CLAS is often used as a "missing mass" spectrometer, where 
all final state particles except one particle are detected. The undetected particle is 
inferred through the overdetermincd kinematics, making use of the good momentum 
and angle resolution. The right panel in figure 3 shows an example of the kinematics 
covered in the reaction ep — * epX. It shows the invariant hadronic mass W versus 
the missing mass Mx- The undetected particles 7r°, rj, and u> are clearly visible as 
bands of constant Mx ■ The correlation of certain final states with specific resonance 
excitations is also clearly seen. 

2.1.3. Experimental Hall C - HMS and SOS 

Hall C houses the high momentum spectrometer (HMS) and the short orbit 
spectrometer (SOS). The HMS reaches a maximum momentum of 7 GeV/c, while 
the SOS is limited to about 1.8 GeV/c. The spectrometer pair has been used to 
measure the 7*NA(1232) and 7* AW* (1535) transition at high Q 2 values. For these 
kinematics the SOS was used as electron spectrometer and the HMS to detect the 
proton. To achieve a large kinematics coverage, the spectrometers have to be moved 
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Figure 3: Left panel: Charged particle identification in CLAS. The reconstructed 
mass/Z(charge number) for positive tracks from a carbon target is shown. Additional 
sensitivity to high-mass particles is obtained by requiring large energy loss in the scintilla- 
tors (shaded histogram). Right panel: Invariant mass versus missing mass for ep — » epX 
at an electron beam energy of 4 GeV. 



in angles, and the spectrometer optics has to be adjusted to accomodate different 
particle momenta. This makes such a two spectrometer setup most useful for study- 
ing meson production at high momentum transfer, or close to threshold. In either 
case, the Lorentz boost guarantees that particles are produced in a relatively nar- 
row cone around the virtual photon, and can be detected in magnetic spectrometers 
with relatively small solid angles. 

2.2. MAMI-B 

The MAMI-B microtron electron accelerator 30 at Mainz in Germany reaches 
a maximum beam energy of 850 MeV. There are experimental areas for electron 
scattering experiments with three focussing magnetic spectrometers with high res- 
olution 31 . A two-spectrometer configuration has been used in cross section and 
polarization asymmetry measurements of n° electroproduction from protons in the 
A(1232) region. 

Another experimental area is equipped for physics with an energy-tagged pho- 
ton bremsstrahlung beam 32 . Experimental setups with BaF-i crystals (TAPS) have 
been employed for measurements of differential cross sections for n° and rj produc- 
tion and for beam asymmetry measurements using a linearly polarized coherent 
bremsstrahlung beam. 

2.3. MIT-Bates 

The Bates 850 MeV linear electron accelerator has been used to study ir° pro- 
duction in the A(1232) region using an out-of-plane spectrometer setup 33 . A set 
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of four independent focussing spectrometers was used to measure various response 
functions, including the beam helicity-dependcnt out-of-plane response function. 
Because of the small solid angles covered by this seteup, a limited range of the 
polar angles in the center of mass frame of the p-K° subsystem could be covered. 
These spectrometers are no longer in use, but data are still being analyzed. 

2.4. Laser backs cattering photon facilities 

Electron storage rings built as light sources for material science studies are often 
used parasitically to produce high energy photons for nuclear physics applications. 
An intense laser beam is directed tangcntially at the electron beam circulating in the 
storage ring producing high energy Compton backscattered photons in an energy 
range dependent upon the wavelength of the laser light. While the photon intensities 
are quite modest, the energy spectrum is peaked at the high energy end providing 
an efficient source of high energy photons for nuclear physics experiments. The 
laser light is easily polarized linearly or circularly. In the Compton backscattering 
process the polarization of the laser light is transferred to the high energy photon 
beam providing a convenient source of polarized photons. 

2.4.1. The Graal Tagged Photon Facility 

The Grenoble Synchrotron Light Source facility is used to generate a laser 
backscattered polarized photon beam of up to 1470 MeV energy for nuclear physics 
applications. A BGO crystal detector is used for the detection of photons 34 cover- 
ing a large portion of 47r. Multi-wire proportional chambers allow charged particle 
tracking. Particle identification is achieved by time-of-flight measurements at for- 
ward angles, and by energy loss measurements at large angles. The large solid angle 
coverage allows the study of reactions with multiple photons in the final state which 
is important for nucleon resonance studies in ir° and rj production 35 . 

2.4.2. The LEGS at Brookhaven National Laboratory 

Brookhaven National Laboratory operates an electron synchrotron as a light 
source with an energy of 2.8 GeV/c. A laser backscattered real photon beam with an 
energy up to 470 MeV is used for nuclear physics experiments 36 . A tagging system 
measures the energy of the Compton-scattered electron from which the photon 
energy is inferred. The photon beam is used with an unpolarized hydrogen or 
nuclear target, and with a polarized HD target 37 . Several arrays of Nal(Tl) crystal 
detectors have been used to measure Compton scattering and tt° production off 
protons in the A(1232) region. 

2.4.3. LEPS at Spring-8 

SPring-8 operates an 8 GeV electron synchrotron near Osaka in Japan. A laser- 
backscattered, energy-tagged polarized photon beam with an energy up to 2.4 GeV 
is produced for nuclear and particle physics applications 38 . The LEPS detector 
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consists of a plastic scintillator to detect charged particles produced in the target, 
an aerogel Cerenkov counter for particle identification, charged-particle tracking 
counters, a large dipole magnet, and a time-of-flight wall for particle identification. 
The LEPS detector has been used for u, and near-threshold <p production, and for 
strange particle production. 

2.5. Electron Stretcher and Accelerator (ELS A). 

The University of Bonn operates a 2.5 GeV electron synchrotron and a stretcher 
ring and post accelerator to obtain a high duty factor beam and an energy of 3 
GeV. Three experimental setups have been used for meson production experiments 
during the past decade. 

2.5.1. The SAPHIR Detector 

SAPHIR is a large acceptance detector with 2n azimuthal coverage. An external 
electron beam was used to generate an energy-tagged real photon beam for exper- 
iments with the SAPHIR detector 39 . At the core of the detector is a large-gap 
dipole magnet. Tracking is provided by a central drift chamber located inside the 
dipole magnet, and additional chambers outside the magnetic field region. Scintil- 
lation counters are used for triggering and to provide time-of-flight information for 
particle identification. 

2.5.2. The Crystal Barrel Detector at ELS A 

The Crystal Barrel (CB-ELSA) detector was originally used at the LEAR pp 
ring at CERN. The detector was recently brought to ELSA for operation in an 
energy-tagged bremsstrahlung photon beam 40 . The detector consists of Csl crystals 
providing nearly full solid angle coverage for neutral particle detection. The main 
focus is the detection of multiple neutral particle final states. 

2.5.3. The Elan Apparatus 

The Elan apparatus has been used for studies of single pion electroproduction 
in the A(1232) region. A focussing magnetic spectrometer detects the scattered 
electrons, and electromagnetic shower detectors measure photons from n° decays. 
Protons and charged pions are detected as well. Charged hadrons are not magnet- 
ically analyzed. This setup has been used for measurements of ir° and n + in the 
A(1232) region. 

3. General Formalism 

The bulk of data from the facilities described in the previous section are from 
experiments with a single meson and baryon in the final state. We therefore only 
present the formulation for such reactions. The generalization of the formulation to 
the cases that the final states are three-body states is straightforward. 

We consider the process N(e, e'M)B illustrated in Fig. 4. The final meson- 
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Figure 4: 
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One-photon exchange mechanism for meson elcctroproduction from a 



baryon states are two-body states, such as ttN, r)N, KA, luN and <j>N . Within the 
Relativistic Quantum Field Theory the Hamiltonian density for describing this 
process can be concisely written as 

H em (x) = eA^x){j»(x) + J»(x)] , (1) 

where is the photon field, 

f{x)=$ e (x) 1 ^ e (x) ( 2 ) 

is the lepton current, and the electromagnetic interactions involving hadrons are 
induced by the hadron current J M . 

With the convention of Bjorken and Drell 41 , the Hamiltonian density Eq.(l) 
leads to 

<kp'\ J dxA' 1 (x)J^(x) | qp >= {2n) A 5 A {p + q - k - p') < hp' \ e^,(q)J^(0) \ qp > , (3) 

where q, p, k, and p' are the momenta for the intial photon, intial nucleon, fi- 
nal meson, and final nucleon, respectively, e^q) is the photon polarization vector. 
Throughout this paper, we will suppress the spin and isospin indices unless they 
are needed for detailed explanations. 
It is convenient to write 



< hp' | e M (g)J'*(0) | qp > = -J- / m ^ / _l_ ^(g)J*(fc,p'; q,p) 



(2wfy E N (p>) ^/2Ejjc) 

-WiT 1 (4) 



E N {p) V2w ' 

The expression for calculating electromagnetic meson production cross sections can 
be expressed in terms of J^{k,p' '; q,p). For evaluating electroproduction cross sec- 
tions, it is common and convenient to choose a coordinate system that the virtual 
photon is in the quantization z— direction, and the angle between the e — e' plane 
and M — B plane is <f>M, as illustrated in Fig. 5. With some straightforward but 
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Figure 5: Kinematics of meson electroproduction reaction. 



lengthy derivations, it is possible to write the differential cross section of N(e, e'M)B 
reaction in the following form 

d 5 (Jh „ d<r h 



dE e id£l e id£l* M 



dn* M 



with 



d<Jh 
dfl M 



= n 



da 



unpol 



da L T' 



+ hy/2e(l - e) — sin<j> M } , 



(5) 



M u M 

where h is the helicity of the incoming electron. The kinematic factors associated 
with the incoming and outgoing electrons are only contained in the following two 
variables 



= {1 



Q 2 

aK H E' 



-tan 
1 



-}- 



2tt 2 Q 2 E e 1 - e ' 



(6) 
(7) 



where Kh = u — Q 2 /2itin is the virtual photon flux, a — 1/137 is the electro- 
magnetic coupling constant, = (oj, q) is the momentum of the photon, and 



q | 2 — uj 2 . The incident and outgoing electron energies are related to 



g M by 



= E P 



EL 



Q 2 = AE e E' e sin 2 ^-. 



(8) 
(9) 



where e is the angle between the incident and outgoing electrons. 

For investigating N* excitations, the kinematics is often characterized by the 
initial 7*iV invariant mass W and Q 2 . For such a choice, the energy transfer is then 
defined by 



W 2 + Q 2 - m 2 N 
2mjv 



(10) 
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The corresponding electron kinematics for conducting experiments with a given 
(Q 2 , W) can then be evaluated by using Eqs.(8)-(9). 

Note that the differential cross sections in the right-hand side of Eq.(5) is defined 
in the center of mass (cm.) frame of the initial j*N and the final MB systems. 
These quantities must be evaluated in terms of the momenta in that cm. frame. For 
the coordinate system chosen as q || z, all momenta needed in calculating da/d£l* M 
must be transformed by a Lorentz boost with (3 = z | q | /(w + mjv)- In terms of 
variables W and Q 2 , a momentum p£ in the considered cm. frame is related to a 
momentum p^ in the laboratory frame by 

= <»> 

P% = P x , (12) 
Pi = P v (13) 

Pi = -^° + ^V. d4) 

Specifically, we have for the virtual photon 

I I mN II 

I Qc I = -j^r | q | , (15) 
um N -Q 2 

^ = — w — • (16) 

It is easy to see that Q 2 = q 2 — uj 2 = q 2 — u 2 

We next present formulae for calculating the cm. differential cross sections in 
the right-hand-side of Eq.(5). The unpolarized cross section is given by 

ddunpoi da T da L da TT , — zda LT , . , 

= + e dnh + e ^7 cos20M + VWT7) ^ C0S ^ M ' (17) 

where and (Tlt are called the transverse, longitudinal, polarization, 

and interference cross sections. These four cross sections and the ddLT'/dQ.* in 



Eq.(5) can be written as 

dap 



—M,3(k c ,p c ;q c ,p c ), (18) 



dn* M q, 

where qj = (W 2 — m 2 N )/{2rriN) — Kh is the effective photon cm. momentum, j3 = 
T, L, TT, LT and LT', and the cm. momenta k c , p' c , q c , and p c can be calculated 
from the corresponding momenta in the laboratory frame by using Eqs.(ll)-(14). 
Obviously qj —\ q c | at the photon point Q 2 = 0. 

The meson production dynamics is contained in Mp of Eq.(18). They are calcu- 
lated from various combinations of current matrix elements evaluated on the 4*m = 
plane (see Fig. 5): 

M T (k c ,p' c ;q c ,p c ) = — ^ [| J x (k c ,p' c ;q c ,p c ) \ 2 + \ J y (k c ,p' c ;q c ,p c ) | 2 ]</, M= o 

spins 
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M L {k c ,p' c ;q c ,p c ) 
M T T(kc,p' c ;qc,Pc) 
M LT (k c ,p' c ;q c ,p c ) 
M LT ,(k c ,p' c ;q c ,p c ) 

with 



^ Yl %tl JZ ( fc c,p' c ;<7c,p c ) | 2 W=o 

spns 

J [I JX ( k c,Pc'iqc,Pc) I 2 - I J y (k C iP f c ;q C iPc) \ 2 Um=o 

spins 

F I Q 2 

spins 

F I Q 2 

V ~ Im { JZ ( k c,P'cl1c,Pc)J X *{kc,Pc;<lc,Pc)}<l>M=0 



2 — v oj 

spins 



(19) 



m N m N E M (k c )E B (p' c ) 



(2Tr) 2 2E M (k c )E B (p' c )E N (p c ) 



2W 



(20) 



where E a (p) — y/p 2 + m 2 with m a denoting the mass of particle a. 

The differential cross sections of N(e,e'M)B are often expressed in terms of 
response functions 42 R a which are related to the differential cross sections of Eq.(18) 

by 



Ri 





1 kc | 




q2 


d&TT 


|k c | 


dn* M 


ql 


da L 


1 k c | 


dn* M 


q2 


da L T 


1 k c 1 


HQ* 


ql 


da L T> 


|k c | 


dn* M 


7 

q2 



Rtt 

Q 2 



Rl 



-R 



LT 



-^jRhT'- 



(21) 

The above formulation can be readily used to calculate various polarization ob- 
servables with a polarized initial nucleon. For observables with a polarized recoiled 
final baryon, the situation is more complicated. They have been explicitly derived 
for pseudo-scalar meson production 43 > 44 . Formulations for analyzing spin observ- 
ables of vector meson production were developed in Ref. 45 . 

We also note that the unpolarized photoproduction cross section is given by 
daT/dVl* M evaluated at Q 2 = and qj — >| q c |. For polarized photons, one needs to 
choose an appropriate combination of J x and J v . For instance, da^/dQ* M (d<j\\/ dVt* M )ior 
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the photon polarization normal (parallel) to the hadron plane is calculated from 
keeping only J y (J x ) contribution and multipling the resulting cross section by a 
factor of 2. The photon asymmetry is defined as 




&± + (T\\ 



(22) 



Calculations of other photoproduction polarization observables are given, for exam- 
ple, in the appendix C of Ref. 66 . 

We next present formulae which are often used in analyzing the production of 
pseudo-scalar mesons, such as MB = nN,r]N, KY . The Lorentz invariance and 
gauge invariance allow us to write the hadron current matrix elements as 



where u(jp) is the Dirac spinor, Ai(s,t,u) are Lorentz invariant functions, and Mi 
are independent invariances formed from 7^, 75, and momenta variables. The 
expressions for Mi are irrelevant to this paper and hence are omitted here. But they 
can be found, for example, on page 5 of Ref. 46 . For ir production, the amplitudes 
defined above can be further classified by isospin quantum numbers. There are 
for the isoscalar photon, and for the isovector the two amplitudes ^t 1 / 2 ) and ^4.( 3 / 2 ) 
for the final wN system with total isospin 1=1/2 and 1 = 3/2 respectively. Each 
invariant amplitude in Eq.(23) can be expanded as 



where r is the isospin Pauli operator, and a is the isospin quantum number asso- 
ciated with the produced pion. Eq.(24) then leads to A^ /2) = + 2A ( f ) and 
A- 3 ^ 2 ' = A^ — A^~\ It is useful to further define proton p^ 1 / 2 ) and neutron 
n^ 1 / 2 ) amplitudes with total isospin 1=1/2 



e ft {q)J' l {k,p';q, P ) = ^ u[{p')A i {s,t,u)M i }u{p) 



(23) 



i=l,6 



^ = ^4 ( ~ ) [^ 3 ] + 4 (+) <5 Q , 3 + 4 



(24) 




4 (0) , I 4(1/2) 

4(0) 1 4(1/2) 

1 3 4 ' 



(25) 



Then the amplitudes for four physical processes can be written as 



A l (j*n 



Ai(7*n 






mr ) 



V2[ P 4 (1/2) - 14 (3/2) ] 

4(1/2) 2 (3/2 ) 

.(1/2) 2 (3/2 ) 
n^i-i 3 * 



(26) 
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The above invariant functions A\ ' are the starting point for developing dispersion 
relation approach which will be given in section 4.7. The isospin relations Eqs.(24)- 
(26) are valid for all of the amplitudes we are going to discuss. However, the 
isospin quantum numbers as well as spin quantum numbers will be suppressed in 
the remainder of this article. 

For investigating nuclcon resonances, it is useful to have a formulation expressing 
the meson production cross sections in terms of multipolc amplitudes. If the final 
hadron state consists of only a pseudo-scalar and a spin 1 /2 baryon, such as wN, KY 
and rjN states, such a formulation has been well developed. This is accomplished 
by casting Eq.(23) into the Chew, Goldbcrger, Low, and Nambu (CGLN) 47 form 
defined in the cm. frame of the final meson-baryon system 

e^J^(k c ,p' c ;q c ,p c ) = ]T F^s, t, Q 2 )u(p c ')OMPc) , (27) 

i=l,6 

where Fi(s,t,u) are the Lorentz invariant CGLN amplitudes and Oj are operators 
defined in the baryon spin space 



with 



0\ = ia ■ b , 


(28) 


O2 = a ■ k c cr • (q x b) , 


(29) 


3 = ia ■ q c k c • b , 


(30) 


O4 = ia ■ k c k c • b , 


(31) 


5 = -ia ■ k c b , 


(32) 


O e = -ia ■ q c 6o , 


(33) 


= e"(«c) - f^j^ 

1 1 


(34) 



Obviously we have b • q c = 0. The CGLN amplitudes Fi(s,t, Q 2 ) can be expanded 
in terms of multipole amplitudes characterized by the angular momentum quantum 
numbers of the initial j*N and the final MB systems. The relations are found to 
be 

F i = £[^+1 0*)^+ + P'l-MEt- + P' l+l {x)Mi+ + (£ + l)PU{x)Me4Zh) 
F 2 = Yy< i+1 )^)M t+ +lP' l {x)M t \ (36) 
F z = ^KjW^ + F/^W^- - + 1 (37) 

p s = ^[-{£+l)Pi{x)S t+ +tPi{x)S t -] (39) 
F 6 = Y,W+ 1 ) P 'i+^ x ) S i+- tP 'i-i{x)St-] (40) 
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In the above equations, the multipolc amplitudes Ei±, Mg± and St± are functions 
of W and Q 2 only. They describe the transitions which can be classified according 
to the character of the photon, transverse or scalar(or longitudinal), and the total 
angular momentum J = £±1/2 of the final state. In addition, the transverse photon 
states can either be electric with parity (— l) Ll , or magnetic, with parity (— 1) L "- +1 , 
where L 7 is the orbital angular momentum of the j*N system. In Table 1 1, wc 
list how each multipolc amplitude with J < 3/2 is related to the initial L 1 and 
final (£, J) angular momentum quantum numbers. The longitudinal multipoles are 
related to the scalar multipoles by Li± = (oj/ | q \)Sg±. 



Table 1: Angular momentum quantum numbers associated with 7* TV — > irN mul- 
tipolc amplitudes. See text for the explanations. 



£ 


J 




Notation 





1/2 


1 


Eo+ 


1 


3/2 


2 


E\+ 


1 


1/2 


1 


Mi_ 


1 


3/2 


1 


Mi+ 





1/2 


1 


Sq+ 


1 


1/2 





Si- 


1 


3/2 


2 


Si+ 



We now note that the matrix elements J % with i = x,y, z for evaluating Eq.(19) 
can be obtained from Eq.(27) by setting e M = (0, x),e^ — (0, y), = (0,5), respec- 
tively. By further using the relations Eqs.(35)-(40), the differential cross sections 
Eq.(5) or Eq.(18) can then be expressed in terms of multipole amplitudes. For 
example, Eq.(5) can lead to the total inclusive cross section 

da , Q 2 . . 

T i a T + t^Io^l] , (41) 



dE' e dQ' e I q | 2 

where 

<?T = -^Y J ^+l)[\M i+ \ 2 + \E (l+1) _\ 2 +f 2 {£+l)[\M t _\ 2 + \E (l _ 1)+ \ 2 ] 

**• = ^tE^+^i^+d-I'+^i^-d+i 2 ]- ( 42 ) 

1 q ° 1 i 

With the above formulation, we then turn to describe various theoretical models 
for analyzing electromagnetic meson production reactions. 



4. Theoretical Models 

The development of theoretical models for investigating electromagnetic pion 
production reactions began in 1950's with the pioneering work by Chew, Goldberger, 
Low, and Nambu(CGLN) 47 . In the subsequent years, their dispersion- relation ap- 
proach was the basis of many analyses 49 of pion production data in the A exci- 
tation region. This approach has been revived 50 ' 51 recently and extended 52 ' 53 to 
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also analyze r\ production. For investigating the data at higher energies where 
the production of two pions and other mesons(77 and K, u), and <f>) could arise , 
the isobar models 48 were developed to extract the parameters of higher mass nu- 
cleon resonances. During the years around 1980, the K-matrix effective Lagrangian 
models 54 ' 55 were developed to study the A excitation. The K-matrix method and 
isobar parameterization have been used subsequently to develop tools for perform- 
ing amplitude analyses of the data and determining the resonance parameters. Ex- 
amples are the very useful dial-in codes SAID 56 and MAID 57 . Progress has also 
been made in extracting resonance parameters using the multi-channel K-matrix 
method 58 ' 59 ' 60 and the unitary coupled-channel isobar model 61 ' 62 ' 63 . 

In recent years, a rather different theortical point of view has been taken to 



develop dynamical mo dels 64 ' 65 ' 66 ' 67 ' 68 ' 69 ' 70 ' 71 ' 72 ' 73 ' 74 ' 75 ' 76 ' 77 ' 78 of meson production 



reactions. These models account for the off-shell scattering effects and can therefore 
provide a much more direct way to interpret the resonance parameters in terms of 
the existing hadron structure models. So far, the dynamical reaction model has 
been able to interprete the resonance parameters, in particular the A resonance, in 
terms of constituent quark models. Its connection with the results from quenched 
and unquenched Lattice QCD calculations remains to be established. 

In the first part of this section, we will give a general derivation of most of the 
exisiting models in order to clarify their differences. We then give some detailed 
formula for the dynamical model which are needed for discussing the results in 
section 5. The analyses based on the dispersion relation approach will be described 
at the end of this section. 

4.1. Hamiltonian Formulation 

Most of the existing models for analyzing the data of electromagnetic meson 
production reactions can be schematically derived from a Hamiltonian formulation 
of the problem. The starting point of our derivation is to assume that the meson- 
baryon (MB) reactions can be described by a Hamiltonian of the following form 



Here v bg is the non- resonant (background) term due to the mechanisms such as 
the tree-diagram mechanisms illustrated in Fig. 6(a)-(d), and v R describes the N* 
excitation in Fig. 6(e). Schematically, the resonant term can be written as 



H = H a + V 



(43) 



where H is the free Hamiltonian and 



V = v bg + v 



(44) 




(45) 



where Ti defines the decay of the z-th N* state into meson-baryon states, and M° 
is a mass parameter related to the resonance position. 
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Figure 6: Tree diagrams for meson-baryon interactions. TV* is a nucleon resonance. 

The next step is to define a channel space spanned by the considered meson- 
baryon (MB) channels: jN, ttN, r/N, 7rA, pN aN 7 ■■. The S-matrix of the meson- 
baryon reaction is defined by 

S(E) a>b = 5 a , b - 2mS(E - H )T a . b (E) , (46) 

where (a, b) denote MB channels, and the scattering T-matrix is defined by the 
following coupled-channel equation 

T a , b (E) = V a . b + J2 V a ,c9c(E)T c , b (E) . (47) 

c 

Here the meson-baryon propagator of channel c is 

g c (E) =< c | g(E) \ c > 

with 

9(E) 



E-Ho+ie 
Pi 



= g'(E)-iTrS(E-H ), (48) 
where 

9 P (E) = (49) 

Here P denotes taking the principal-value part of any integration over the propa- 
gator. We can also define K-matrix as 

K a . b (E) = V a , b + V a ^(E)K c . b (E) . (50) 

C 

Eqs.(47)-(50) then define the following relation between the K-matrix and T-matrix 
T a , b (E) = K a . b (E) -J2T a , c (E)[in6(E - H )] c K C}b (E) . (51) 
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By using the two potential formulation 94 , one can cast Eq.(47) into the following 
form 

T a . b (E) = t%(E)+t* b (E) (52) 

with 

tU E ) = E %, a {E)[G{E)\ i<j T N ., b {E). (53) 

N*,N* 

The first term of Eq.(52) is determined only by the non- resonant interaction 

t^E)=v% + Y,< g ME)t%{E). (54) 

C 

The resonant amplitude Eq.(53) is determined by the dressed vertex 

Tn-AE) = r N ,, a + J2^N^ b gb(E)tl 9 a (E), (55) 

b 

and the dressed propagator 

[G(E)-\ 3 {E) = (E- M^Sij - Hij(E) . (56) 
Here M%, is the bare mass of the resonance state N* , and the self-energy is 

S W (E) = r]y* ,a9a (E)T N ; ,„ (E) . (57) 
a 

Note that the meson-baryon propagator g a (E) for channels including an unstable 
particle, such as 7rA, pN and aN, must be modified to include a width due to their 
decay into ttttN channel. In the Hamiltonian formulation, this amounts to the 
following replacement 

^)^<a\ E _ Ho 1 _^ v[E) \a>, (58) 

where the energy shift is 

Here Ty describes the decay of p, a or A in the quasi-particle channels. 

Eq.(47), Eqs.(52)-(59), and Eq.(51) are the starting points of our derivations. 
^From now on, we consider the formulation in the partial-wave representation. The 
channel labels, (a, b, c), will also include the usual angular momentum and isospin 
quantum numbers. 
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4.2. Tree-diagram models 

The tree-diagram models are based on the simplification that T ~ V = v b9 +v R . 
The resonant effect is included by rnodifmg the mass parameter of v R , defined in 
Eq.(45), to include a width, such as M° = M i - §rf *(£). Eq.(47) is then simplifed 
into 

where v bg is calculated from the tree-diagrams (Fig. 6(a)- (d)) of a chosen Lagrangian, 
and r* ot is the total decay width of the i-ih N*. 

In recent years, the tree-diagram models have been applied mainly to investi- 
gate the photoproduction and electroproduction of K mesons 79 ' 80 ' 81 ' 82 ' 83 ' 84 , vec- 
tor mesons 85 ' 86 ' 87 (u, <p) and two pions 90 . At high energies, the t-channel am- 
plitudes(Fig. 6(b)-(c)) are replaced by the Regge parameterization in some tree- 
diagram models 88,89 . The validity of using the tree-diagram models to investigate 
nucleon resonances is obviously very questionable, as discussed in a study of u> 
photoproduction 87 and kaon photoproduction 74 . 



4.3. Unitary Isobar Models (UIM) 



4.3.1. MAID: 

The Unitary Isobar Model developed 57 by the Mainz group is based on the on- 
shell relation Eq.(51). By including only one hadron channel, ttN ( or r]N ), Eq.(51) 
leads to 

TttN.iN = TT 7 !? KttN^N 

= e lS7,N cosSknKttN^n ■ (61) 

Here we have used the relation K^n.ttN = —tan5 W N with S^n being the pion-nucleon 
scattering phase shift. By further assuming that K = V = v bg + v R , one can cast 
the above equation into the following form 



T vN „ N { UIM) = e 5 ^cos5Mv%, lN ] + J2 T S 7 jv(^) • (62) 

N* 



Clearly, the non-resonant multi-channel effects, such as — > (pN,irA) — > ttN, 
which could be important in the second and third resonance regions are neglected 
in MAID. In addition, they calculate the non-resonant amplitude v v 9 N N using an 
energy-dependent mixture of PV and PS (pseudo-scalar) ttNN coupling 



A 2 2 

j- _ A m T PV q TPS / fi o\ 

— . 2 T 2 n TrNN ' A 2 _i_ 2^tNN > \ 00 ) 



20 



where qo is the on-shell photon momentum. With cutoff A m = 450 MeV, one then 
gets PV coupling at low energies and PS coupling at high energies. 

For resonant terms in Eq.(62), MAID uses the following Walker's parameterization 48 

tS„ n (E) = fUE) Mf ^2^^ W E ) Ai > ( 64 ) 

where f^ N (E) and f* N (E) are the form factors describing the decays of N*, T tot 
is the total decay width, A 1 is the -fN — > A^* excitation strength. The phase $ is 
determined by the unitary condition and the assumption that the phase ip of the 
total amplitude is related to irN phase shift S^n and inelasicity tj^n by 

= l- V ME)cos25ME) 
K ' L ^ N {E)sin25„ N {E) J 1 ; 

4.3.2. JLab/Yeveran UIM: 

The Jlab/ Yerevan UIM 53 is similar to MAID. But it implements the Regge 
parameterization in calculating the amplitudes at high energies. It also uses a 
different procedure to unitarize the amplitudes. 



Both MAID and JLab/Yeveran UIM have been applied extensively to analyze 
the data of ir and 77 production reactions, as will be discussed in section 5. Very 
useful new information on A^* have been extracted. 

4.4. Multi-channel K-matrix models 

AAA. SAID: 

The model employed in SAID 56 is based on the on-shell relation Eq.(51) with 
three channels: jN, irN, and ttA which represents all other open channels. The 
solution of the resulting 3x3 matrix equation can be written as 

T~/n,kn(SAID) = A/(l + iTkn^n) + A R T vN}7vN , (66) 

where 

A ts K~{N,irAK-!TN,TrN , an \ 

m = A 7 jv,7riV 77 , [pi) 

R--kN,-kA 

An = (68) 
In actual analyses, they simply parameterize Aj and Ar as 

M 

A I = [v%, v N] + ^2PnZQl a +nC), ( 69 ) 
n=0 
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A R = !^(f^f>„(^ r , (70) 

where ko and go are the on-shell momenta for pion and photon respectively, z = 
^/fcg + 4m^/&o, Ql{z) is the legendre polynomial of second kind, E n = E 1 — m^{l-\- 
■m^j (2m7v)), and p n and p n are free parameters. SAID calculates v 9 N ^ N of Eq.(69) 
from the standard PS Born term and p and u> exchanges. The empirical ttN ampli- 
tude T^jv^at needed to evaluate Eq.(66) is also available in SAID. 

Once the parameters p n and p n in Eqs.(69)-(70) are determined, the N* param- 
eters are then extracted by fitting the resulting amplitude T^m^n at energies near 
the resonance position to a Breit-Wigner parameterization(similar to Eq.(64)). Very 
extensive data of pion photoproduction have been analyzed by SAID. The extension 
of SAID to also analyze pion electroproduction data is being pursued. 

4.4.2. Giessen Model 

The coupled-channel model developed by the Giessen group 60 can be obtained 
from Eq.(51) by taking the appoximation K = V; namely, neglecting all multiple- 
scattering effects included in Eq.(50) for K-matrix. This leads to a matrix equation 
involving only the on-shell matrix elements of V 

T a , b (Giessen) = ^[(1 + iViE))- 1 ^^) . (71) 

c 

The interaction V = v bg + v R is evaluated from tree-diagrams of various effective 
Lagrangians. The form factors, coupling constants, and resonance parameters are 
adjusted to fit both the ttN and 7^ reaction data. They include up to 5 channels in 
some fits, and have identified several new N* states. But further confirmations are 
needed to establish their findings conclusively, as will be discussed later in section 
5.6. 



4.4.3. KSU Model 

The Kent State University (KSU) model 58 can be derived by noting that the 
non-resonant amplitude t b9 , defined by a hermitian v bg in Eq.(54), define a S-matrix 
with the following properties 

S b a 9 b (E) = 6 a , b -2ni6(E-H )t b a 9 b (E) (72) 

= ^t )T (E)J c ^(E), (73) 

C 

where the non-resonant scattering operator is 

uitXm = 5 a , c +g a (E)t%(E). (74) 

With some derivations, the S-matrix Eq.(46) and the scattering T-matrix defined 
by Eqs.(52)-(57) can then be cast into following form 

S a , b (E) =J2^ T (E)RcAE)^(E) , (75) 

c,d 
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with 

R c>d (E)=S c , d + 2iT R d (E). (76) 

(77) 

Here we have defined 

T c R d (E) =Y,^N : , c m G ( E ^ T ^^E) . (78) 

The above set of equations is identical to that used in the KSU model of Ref. 58 . 
In practice, the KSU model fits the data by parameterizing T R as a Breit-Wigner 
resonant form T R ~ xV/2/(E - M - iT/2) and setting = B = B X B 2 ■ ■ ■ B n , 
where Bi = exp(iXAi) is a unitary matrix. 

The KSU model has been applied to ttN reactions, including pion photoproduc- 
tion. It is now being extended to investigate KN reactions. 



4.5. The CMB Model 

A unitary multi-channel isobar model with analyticity was developed 61 in 1970's 
by the Carnegie-Mellon Bcrkeley(CMB) collaboration for analyzing the 7riV data. 
The CMB model can be derived by assuming that the non-resonant potential v bg 
is also of the separable form of v R of Eq.(45) 

b rl a r L , b r^ a T H . b 

a ' b E — Ml E — Mn [ ' 

The resulting coupled-channel equations are identical to Eqs.(52)-(59), except that 
t b a 9 b = and the sum over N* is now extended to include these two distance poles 
L and H . 

By changing the integration variables and adding a substraction term, Eq.(57) 
for the self-energy can leads to CMB's dispersion relations 

^A 3 ) = ^li,c<&c{s)lj,c, (80) 

C 

Re[Z c (s)] = Re[* c (so)]+ S — T , / m[ f e( / )] . ds', (81) 

77 Js th ( s - s )( s - s 0) 

where 7i jC is a coupling constant defining the decay of N* into channel c. Thus 
CMB model is analytic in structure which marks its difference with all K-matrix 
models described above. 

The CMB model has been revived in recent years by the Zagreb group 62 and a 
Pittsburgh- ANL collaboration 63 to extract the N* parameters from fitting the re- 
cent empirical wN and 7./V reaction amplitudes. The resulting N* parameters have 
very significant differences with what are listed by PDG in some partial waves. In 
particular, several important issues concerning the extraction of the N* parameters 
in Sn channel have been analyzed in detail. 
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4.6. Dynamical Models 



A. In the A region 

Keeping only one resonance N* — A and two channels a,b — irN, jN ,Eqs.(52)- 
(57) are reduced to what were developed in the Sato-Lee (SL) model 68 ' 70 . 




Quark Model Prediction 



+ <=sa — i — *- 

N, A 

1 



* 

pion cloud effect 

Figure 7: Graphic representation of the dressed -fN — * A vertex. 
Explicitly, we have 

b ^ A7rN (E)T AtWN (E) 
T v n !wN {E) = t n 9 N ^ N (E) + E _ M ^_ z A ( E j > ( 82 ) 

rl t7N (E)T A ,„ N (E) 

•a 



T lN ,ME) - t 7 9 w (E)+ e _ m o a _^ a{e) . (83) 



with 



r A ,j N (E) = r AtjN + T AtnN G nN (E)tl 9 N ^ N (E), (84) 

Fa,ttn(E) = T At7rN + T Ai7rN G 7r N(E)t b 7V 9 N ^ N (E) , (85) 

t a N, lN {E) = v% !7N +t%^ N (E)GME)v%^ N , (86) 

t%,„ N (E) = v% tVN + v% tVN GME)t%^ N (E), (87) 



and 



X A (E) - T AvN (E)GME)T A ,ME)- (88) 



The above equations clearly indicate how the non-resonant interaction modify the 
resonant amplitude. Specifically, Eq.(84) for the dressed A — > 7^ within the SL 
model is illustrated in Fig. 7. 

Alternatively, we can also cast Eq.(47) in the A region as 

T jNt7rN (E) = tf N ^ N (E) + t^ N7lN (E) (89) 
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with 

t^ N (E) = v% tVN + v b 7 % !wN GME)T*N,ME) (90) 
t*N,nN( E ) = v iN,itN + Vj Nt7rN G nN (E)T 7vNt7rN (E) (91) 

The above equations are used by the Dubna-Mainz- Taiwan (DMT) model 71,73 ex- 
cept that they depart from a consistent Hamiltonian formulation and replace the 
term t^ N N by the Walker's parameterization 48 

^n.^n(E) = Un(E) M 2 ^£ 2 MA ^ MArtot fjN(E) A lN . (92) 

Other differences between the SL Model and the DMT model are in the employed 
ttN potential and how the non-resonant jN — > ttN amplitudes are regularized. In 
the DMT model, the non-resonant 7^ — > ttN amplitudes are calculated by using 
MAID's mixture Eq.(63) of PS and PV couplings, while their irN potential is from a 
model 95 using PV coupling. In the SL model, the standard PV coupling is used in a 
consistent derivation of both the ttN potential and 7iV — > 7riV transition interaction 
using a unitary transformation method. 

We now turn to giving relevant formula which are needed for our discussions 
in section 5.1 on the A resonance. The A excitation is parameterized in terms of 
Rarita-Schwinger field. In the A rest frame where toa = qo + EN(q), the resulting 
jN — > A vertex function can be written in the following more transparent form 



, . e j E N (q) + m N 1 3(m A + m N ) 

<A|1 7A ^a| 9 > - (27r )3/ 2 y 2EN{q) ^Am N (E N (q) + m N ) 13 



[iG M {q z )S x q ■ e + G E (q')(S ■ ecr ■ q + S ■ qa ■ e) 
+^£kls . qa . qe ], (93) 

TOA 



where e = 47r/137, q = (uj,q) is the photon four- momentum, and e = (eo, e) is the 
photon polarization vector. The transition operators S and T are defined by the 
reduced matrix element <A||S ? ||7V>=<A||T||iV>=2in Edmonds' convention 96 . 
By using Eq.(93) and the standard definitions 97 ' 98 of the multipolc amplitudes, 
it is straightforward to evaluate the magnetic Ml, electric E2 and Coulomb C2 
amplitudes of the jN — > A transition. We find 70 that 

G M {q 2 ) = ^[I>^a]mi, (94) 

G E (q 2 ) = ^[T iN ^a]e2, (95) 

Gc(q 2 ) = |^[r 7J v^A]c2, (96) 

with 

/toa \q\ 1 



N = 



2m n V mN [1 — q 2 / (m^ + m^) 2 ] 1 / 2 
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At q 2 = 0, the above relations agree with that given in Appendix A of Ref. 66 . 
Equations (93)-(96) can also be used to relate the dressed vertex F 7 jv^A, defined 
by Eq.(84), to the corresponding dressed form factors: 

G * M {q 2 ) = ■^[T 7 jv->a]mi ) 

G* E {q 2 ) = "^-[f 7 AT^A]£;2, 
G c(9 2 ) = |-Tj-^[f 7 Ar^A]c2- 

At the q 2 — photon point, we will also compare our results with the helicity 
amplitudes defined by PDG 24 . They are related to the multipole amplitudes defined 
above by 

V3 

^3/2 = — [[r 7 Af^A]_E2 — [r 7 jV->A]Ml], (97) 

■A-i/2 = — 2 [3[f 7 jv^A]E2 + [T 7 jv->a]mi]- (98) 

At the A resonance position E = Mr = 1236 MeV, the 7riV phase shift in the 
P 33 channel goes through 90 degrees. This leads to a relation, as derived in detail 
in Ref. 68 , that the multipole components of the dressed vertex f 7 ^r^A are related 
to the imaginary(im) parts of the 7./V — > ttN multipole amplitudes in the nN P 33 
channel 



2\ 1 \rK 1 1 8^m A kT A / fl/f 3/2x 



o-M) - ^^.^I^xmO (ioo) 



2^ Gc(9) = 7V [r ^ A]g2= jvV 3m wg x Im ^i), (101) 

where Ta is the A width, fc and g are respectively the momenta of the pion and 
photon in the rest frame of the A. Note that the upper index K in T^^ N in 
Eqs.(99)-(101) means taking only the principal-value integration in evaluating the 
second term of Eq.(84). Details are discussed in Ref. 68 . 

^From the above relations, we obtain a very useful relation that the E2/M1 
ratio Rem and C2/M1 ratio Rsm of the dressed 7^V — > A transition at W = 1232 
MeV can be evaluated directly by using the 7iV — > ttN multipole amplitudes 

[T^a]£2 _ Im(< 2 ) no2 , 
" [r^Al^-im^O' (1 ° 2) 

[f 7 V^ A ]c2 Im(^( 2 ) 
flSM — 7^7? i — TFT- U U <J) 
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Eqs.(99)-(103) can be used in the empirical amplitude analyses to extract the 
form factors and the E2/M1 and C2/M1 ratios of the 7./V — > A transition. The 
extractions of the bare vertices, which can be compared with the predictions from 
most of the constituent quark model calculations, can only be achieved by using 
the dynamical model through Eq.(84). This indicates why an appropriate reaction 
theory is needed in the N* study, as illustrated in Fig.l. 



B. In the Second and third resonance regions 

In these regions, we need to include more than 7riV channel to solve Eq.(47) or 
Eqs.(52)-(59). In addition, these formula must be extended to include explicitly 
the irirN channel, instead of using the quasi two-particle channels 7rA, pN, and 
aN to simulate the irirN continuum. This however is still being developed 76 . Here 
we continue to explain the current investigations in the second and third resonance 
regions within the formulation defined by Eqs.(52)-(59). 

Eqs.(52)-(59) are used in a 2-N* and 3-channel (ttN, r)N, and ttA) study 69 of 
ttN scattering in Sn partial wave, aiming at investigating how the quark-quark 
interaction in the constituent quark model can be determined directly by using the 
reaction data. Eqs.(52)-(59) are also the basis of examining the N* effects 86 and 
one-loop coupled-channel effects 87 on u> meson photoproduction and the coupled- 
channel effects on K photoproduction 74 . 

The coupled-channel study of both ttN scattering and 7 TV — + ttN in Sn channel 
by Chen et al 72 includes irN, rjN , and 77V channels. Their ttN scattering calculation 
is performed by using Eq.(47), which is of course equivalent to Eqs.(52)-(59). In 
their 7 TV — > ttN calculation, they neglect the 7./V — > rjN — > irN coupled-channel 
effect, and follow the procedure of the DMT model to evaluate the resonant term 
in terms of the Walker's parameterization (Eq.(64)). They find that four TV* are 
needed to fit the empirical amplitudes in Sn channel up to W = 2 GeV. 

A coupled-channel calculation based on Eq.(47) has been carried out by Julich 
group 104 for irN scattering. They are able to describe the wN phase shifts up 
to W = 1.9 GeV by including irN, r]N, 7rA, pN and aN channels and 5 TV* 
resonances : P 33 (1232), 511(1535), Sii(1530), Sii(1650) and D 13 (1520). They 
find that the Roper resonance Pn(1440) is completely due to the meson-exchange 
coupled-channel effects. 

A coupled channel calculation based on Eq.(47) for both 7riV scattering and 
7iV — > ttN up to W — 1.5 GeV has been reported by Fuda and Alarbi 75 . They 
include irN, jN, rjN, and irA channels and 4 N* resonances : P 33 (1232), Pn(1440), 
5n(1535), and Z?i 3 (1520). The parameters arc adjusted to fit the empirical multi- 
pole amplitudes in a few low partial waves. 

Much simpler coupled-channel calculations have been performed by using sepa- 
rable interactions. In the model of Gross and Surya 67 , such separable interactions 
are from simplifying the meson-exchange mechanisms in Figs 6.(a)-c) as a contact 
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term like Fig. 6(d). They include only irN and 77V channels and 3 resonances: 
P 33 (1232), Pn(1440) and £>i 3 (1520), and restrict their investigation up to W < 1.5 
GeV. To account for the inelasticities in Pn and D i3 , the N* — ► ttA coupling is 
introduced in these two partial waves. The inelasticities in other partial waves are 
neglected. 

A similar separable simplification is also used in the chiral coupled-channel 
models 77 ' 78 for strange particle production. There the separable interactions are 
directly determined from the leading contact terms of SU(3) effective chiral La- 
grangian and hence only act on s-wave partial waves. They are able to fit the total 
cross section data for various strange particle production reaction channels without 
introducing resonance states. It remains to be seen whether these models can be 
further improved to account for higher partial waves which are definitely needed to 
give an accurate description of the data even at energies near production thresholds. 

4.7. Dispersion-relation approaches 

Historically, the approach based on dispersion relations is defined within the 
S-matrix theory which was introduced as an alternative to the relativistic quantum 
field theory in investigating non-perturbative hadron interactions. This approach 
was first applied to investigate pion photoproduction by Chew, Goldberger, Low, 
and Nambu 47 (CGLN) and electroproduction by Amaldi, Fubini and Furlan 98 . It 
was fully developed 49,105 in the years around 1970 to analyze the data at energies 
near the A resonance. In recent years, it has been revived by Aznauryan 51,52 , and 
by Hanstein, Drcchsel, and Tiator 50 for investigating pion photoproduction and r\ 
photoproduction 53 . 

The dispersion relation approach assumes that the scattering amplitude is uni- 
tary and possesses various established symmetry properties such as Lorentz invari- 
ance and gauge invariance. The dynamics is defined by the assumed analytical 
property and crossing symmetry. For 7r and r\ production, the starting point is the 
fixed t dispersion relation 106 for the invariant amplitudes Ai defined in Eq.(24) 

Re[A{{s,t)\ = A^ ole (s,t) + -P / ds'[— + ^}Im[A{(s',t)} , (104) 

77 Jsthr s s s u 

where A 1 ^ ^ (s,t) is calculated from pseudo-scalar Born term, 1 = 0,+,— denote 
the isospin component, £i = £ 2 = —£3 = —£4 = 1 and e+ = e° = — e~ = 1 
are defined such that the crossing symmetry relation A T k (s,t, u) = ^e 7 A\.{u,t, s) 
is satisfied. With the definitions Eqs.(23) and (27) and the multipole expansion 
defined by Eqs.(35)-(40), the above fixed-t dispersion relation leads to the following 
set of coupled equations relating the real part and imaginary parts of multipole 
amplitudes 

P f°° 

Re[Ml{W)] = Mj; pole {W) + - / dW V K^{W, W')Im[Ml(W)] , (105) 

77 JW thT J, 

where M\ is the multipole amplitude, M^' pole (W) is calculated from pseudo-scalar 
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Born term, and K\ e , contains various kinematic factors. In the recent work of 
Ref. 50 , the procedures of Ref. 105 are used to solve the above equations by using the 
method of Omnes 107 . It assumes that the multipole amplitude can be written as 

M}(W) = exp i *'W — M\{W), (106) 
rei 

where is a real function and ra is some kinematic factor, and hence 

Im[M}(W)} = hl*(W)Mi(W) , (107) 
with hg(w) — sin(^)exp l ^(^ w \ The phase 4>\ is assummed to be 
,i,tx,n ,\ - ■n I AW)cos28 I AW) . 

where <5j and r\\ are the phase and inelasticity of 7riV scattering in the partial- wave 
with quantum numbers (£,I). 

The next approximation is to limit the sum over I' in the right-hand side of 
Eq.(105) to a cutoff l m ax- For investigating production below E 1 = 500 MeV, 
(-max — 1 is taken. Another approximation is needed to handle the integration 
over W in Eq.(105). In Ref. 50 , the integration is cutoff at W = A = 2 GeV such 
that all needed phase can be determined by the empirical ttN phase shifts. The 
neglected contribution from W > 2 GeV is then acounted for by adding vector 
meson exchange terms. Eq.(105) then becomes 

1 f A h^iW^MKW^dW 



mkw) = M i e > po[e (w)+- r 



w thr W'-W-ie 



+ 



-V / dW'Kl'/,' (WW')hl',*(W') + M I / V (W) (109) 

n e > I' J W thr 



The method for solving Eqs.(109) is given in Ref. 50 . With the above procedures, 
the model contains 10 adjustable parameters. Excellent fit to all 7^ — > nN data 
up to E 7 = 500 MeV has been obtained in Ref. 50 . 

The calculation in Ref. 51 follow the same approach with additional simplification 
that the coupling between different multipoles and the contribution from W > A to 
the integration are neglected; setting K^'|, = and M^ ,y = in solving Eq.(109). 
These simplifications are justified in calculating the dominant A excitation ampli- 
tude M-l + ' ; . But it is questionable if they can be applied for calculating weaker 
amplitudes. Thus no attempt was made in Refs. 51 ' 52 to fit the data directly using 
dispersion relations. Rather, the emphasis was in the interpretation of the empiri- 
cal amplitudes Mj 3 / 2 , E^ 2 in terms of rescattering effects and constituent quark 
model prediction. By assuming the multipole expansion is also valid in electropro- 
duction, the Q 2 -dependence of these A excitation amplitudes are then predicted. 
There are questions regarding the validity of multipole expansion at E 1 > 500 MeV 
and large Q 2 108 . 
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The dispersion relation approach is also used in Ref. 52 to analyze the pion pho- 
toproduction and electroproduction data in the second and third resonance region. 
It is assumed that the imaginary parts of the amplitudes in Mjv + < W < 2 
GeV are from the resonant amplitudes parameterized as the Walker's Breit-Wigncr 
form Eq.(64), and in 2.5 GeV < oo from Regge-pole model. The imaginary part of 
the amplitude in 2 GeV < W < 2.5 GeV is obtained by interpolation. The real part 
of the amplitude is then calculated from the dispersion relation described above. 
The empirical amplitudes are then fitted by adjusting the resonant parameters. It 
turns out that the resulting parameters are close to what were determined in the 
single channel K-matrix model described in subsection 4.3.2. 

With appropriate modifications, the dispersion relation approach can be applied 
to investigate the production of other pseudo-scalar mesons. This has been achieved 
in Ref. 53 in analyzing the data of r\ production reactions. 

5. Data and Results of Analyses 

A large volume of data of electromagnetic meson production reactions is needed 
to extract the fundamental physics on resonance transition form factors or discover 
new baryon states. Efforts in this direction in the 1970's and 80's at various labo- 
ratories were hampered by the low duty cycle synchrotrons that were available for 
these studies, and by the use of magnetic spectrometers with relatively small accep- 
tance. For a discussion on these results see the excellent review by F. Foster and 
G. Hughes 109 . The construction of CW electron accelerators, and the advances in 
detector technologies have made it possible to use detector system with nearly An 
solid angle coverage, and the ability to operate at high luminosity. Moreover, the 
detection of multiple photons from tt° or r\ decays with high resolution has become 
feasible with the development of high density crystals with sufficient light output, 
such as BGO, Csl, and PbF 2 . These detectors have become powerful tools in the 
study of baryon spectroscopy and structure. In this section we will highlight the 
data and review the results from the analyses. We will only consider meson produc- 
tions from nucleon targets. An extensive review of meson production from nuclei 
has been published recently no . 

5.1. Single pion production. 

Single pion photoproduction and electroproduction have been the main processes 
in the study of the electromagnetic transition amplitudes of the lower mass nu- 
cleon resonances such as A(1232), Pn(1440), L>i 3 (1520), 511(1535), and Fi 5 (1680). 
Data on pion photoproduction now exist from LEGS 111 and MAMI 112 , and from 
GRAAL 113 ' 35 j including results from measurements of polarized beam asymmetries 
and beam-target double polarization observables 114 > 115 . 

During the past few years high statistics data of pion electroproduction have 
been collected at JLab. The ep — > epir° data cover a large range in invariant mass W 
from threshold to 2.5 GeV, a wide range in momentum-transfer Q 2 = 0.1 — 6 GeV 2 , 
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and the full range in azimthal and polar angles in the pir° center-of-mass. 

In the past there have been very limited data on mr + , mostly at forward center 
of mass angles 120,121 ' 122 , some at backward angles 123 . Even fewer data exist in 
7r~ production from deuterium 124 . This has limited our ability to extract reliable 
resonance transition amplitudes in the high mass region where many isospin i 
states exist, which couple more strongly to mr + than to pir . New mr + data from 
CLAS 125 have nearly full angular coverage and span the range W = 1.1 - 1.6 GeV 
and Q 2 = 0.3 — 0.6 GeV 2 . They vastly increase the covered kinematics with high 
statistics, and will be extended to W < 2.5 GeV 2 and Q 2 = 0.1-6 GeV 2 . 

For the first time there are also significant amounts of polarized beam asym- 
metries data and data on the beam helicity response function ctlt' available, both 
for pn° !26, 127,128, 129 an( j £ Qr nn + 130^ -phe most complete data sets will come from 
JLab for both the ep — > epir , and the ep — ► enir + channels. Some response func- 
tions at few low Q 2 have also been measured MIT-Bates and Mainz. In particular, 
the experiments using the OOPS of MIT-Bates yield rather precise data of Att 
and Atl at Q 2 = 0.126 (GeV/c) 2 covering a limited angular range. In table 2, we 
summarize these new data. 



Table 2. Summary of the single pion electro-production data 



Reaction 


Observable 


W range 


Q A range 


Lab/ experiment 






(GeV) 


(GeV 2 ) 




ep — » epir° 


dcr/dQ 


< 1.8 


0.4 - 1.8 


JLab-CLAS 126 




dcr/dQ 


A(1232) 


0.1 - 0.9 


ELSA-Elan 131 




dcr/dQ 


A(1232) 


2.8, 4.0 


JLab-Hall C 117 




dcr/dQ 


< 2.5 GeV 


2 - 6 


JLab-CLAS 132 




dcr/dQ 


A(1232) 


7.5 


JLab-CLAS 133 




dcr/dQ 


< 2.0 


1.0 


JLab-Hall A 118 


ep — » en7r+ 


dcr/dQ 


< 1.6 


0.3 - 0.65 


JLab-CLAS 125 




dcr/dQ 


A(1232) 


0.1 - 0.9 


ELSA-Elan 131 




dcr/dQ 


< 2.5 GeV 


2 - 6 


JLab-CLAS 132 


ep — > epir° 


A e , Att, Atl 


A(1232) 


0.2 


MAMI-A1 128 




cr LT i 


A(1232) 


0.3 - 0.65 


JLab-CLAS 127 




A e 


A(1232) 


0.126 


Bates-OOPS 129 


ep — > en7r+ 


cr LT ' 


< 1.6 


0.3 - 0.65 


JLab-CLAS 130 


ep — > epir 


pol. resp. fct. 


A(1232) 


1.0 


JLab-Hall A 134 


ep — > epir° 


Ap , A e p 


A(1232) 


0.5 - 1.5 


JLab-CLAS 136 


ep — » enn + 


A e p 


< 1.85 


0.4, 0.65, 1.1 


JLab-CLAS 135 



One of the main outcomes from the analyses of these single pion production data 
is a more detailed understanding of the A(1232) resonance. The focus has been on 
the determination of the magnetic Ml, electric E2, and Coulomb C2 form factors 
of the 7^V — > A transition. This development will be discussed in detail in this 
subsection. The single pion production in the second and third resonance regions 
will be covered mainly in section 5.3 where some N* parameters extracted from a 
combined analysis including the data of r] production will be discussed. 



31 



5.1.1. Pion photoproduction 

The high statistics of the photon asymmetry data is essential in determining the 
small E±+ amplitude of the 7-/V — > irN reaction, which determines the electric E2 
strength of the 7^ — + A transition through Eq.(lOO). Fig. 8 shows the comparison 
of the results from the Sato-Lee (SL) model 68 and the -yp — ► pn° data from Mainz 
and LEGS. When the E\ + amplitude is turned off in the SL model, the predicted 
photon asymmetries (dotted curves) deviate from the data. By performing the 
amplitude analyses of these new data by several groups, we now have a world 
averaged value of the R EM ratio, defined by Eq.(102), R EM = (-2.38 ± 0.27)% 116 
at photon point. The magnetic Ml transition strength, defined by Eq.(99), has also 
been determined as G* M (0) = 3.18 ± 0.04. 

30 
20 
10 


30 
20 
10 


30 
20 
10 



200 300 400 200 300 400 

Figure 8: The differential cross section(<ier/<if2) and photon asymmetry(E) of the 
p(j, n°)p reaction, calculated from the Sato-Lee model, are compared with the data 
of Mainz 112 . The photon asymmetry data of LEGS 111 agree with those shown here 
and hence are not displayed. The dashed curves are obtained from setting Rem = 0. 
The dashed and solid curves for da/dSl are indistingushable. 

For the dynamical models, it is possible to also get the bare transition strengths 
Gm(0) and G_e(0) which are obtained by separating the pion cloud effects from 
the full(dressed) transition strengths, as defined by Eq.(84) and illustrated in Fig. 7. 
In Table 3, we show the importance of the pion cloud. We see that the helicity 
amplitude ^3/2 extracted from three different analyses are very close to each other 
and are about 40 % larger than the bare strengths extracted within the dynamical 
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Table 3: Helicity amplitude A 3 / 2 and E2/M1 ratio Rem for the — ► A transition at 
Q 2 = photon point. A 3 / 2 is in unit of W^^GeV^ 1 ^ 2 and Hem in %. The references are: 
a(68), b(71), c(55), d(50), e(101), f(19). * This value differs from the value 226 listed in 
Ref.69, because of a kinematic factor is not included correctly. 





A 3/2 
Dressed 


Bare 


Rem 
Dressed 


Bare 


Rcfs. 


Dynamical Model 


-258* 


-153 


-2.7 


-1.3 


a 




-256 


-136 


-2.4 


0.25 


b 


K-Matrix 


-255 




-2.1 




c 


Dispersion 


-252 




-2.5 




d 


Quark Model 




-186 




~ 


e 






-157 




~ 


f 



models of Rcfs. 68,71 . We now note that these bare values are within the ranges 
predicted by two constituent quark models. This suggests that the bare parameters 
of the dynamical model are more likely to be identified with the current hadron 
structure calculations. In Table 3 we also see that the differences between dressed 
and bare values of Rem are even larger. The bare values from two dynamical 
model analyses are quite different, indicating some significant differences in their 
formulations as discussed in section 4. 

5.1.2. Pion electroproduction 

As can be seen in Table 2, pion electroproduction data are now very extensive 
and of high quality. In Figs. 9, 10, and 11, we show some sample data from CLAS 
at JLab. As an example for a spectrum with high statistics data on tt° electro- 
production at a fixed backward angle of = 170° we show in Fig. 12 response 
functions recently obtained from JLab Hall A 118 . 

In Fig. 9 and Fig. 11, the predictions from the SL, MAID, and DMT models are 
also displayed to illustrate the status of current reaction models. The analyses of 
these new data in the past few years have led to rather accurate determinations of 
the 7-/V — > A transition form factors. We now discuss this advance in more detail. 

5.1.3. The jN -> A(1232) transition form factors 

With the fairly extensive coverage over angles and energies, the data from JLab 
have allowed nearly model-independent determinations of 7A — ► A(1232) form 
factors. Theses analyses by the CLAS collaboration are based on the following 
considerations. At the A peak, the dominant amplitude is M\ + and the small E\ + 
and S\+ can become accessible through their interference with the dominant M\ + 
amplitude. One thus can start the analysis by using a truncation, in which only 
terms involving M\ + are retained. With the partial-wave decomposition defined by 
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Figure 10: ep —* eri7r data from CLAS. Separated response functions or + eo"L, cttt, 
and o\lt are shown. The curve represents a fit to the data using the JLab/ Yeveran unitary 
isobar model . 
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Figure 11: CLAS data on ult 1 °{p( e , e'n + )n reaction in the A(1232) region are com- 
pared with predictions from SL(solid), MAID(dashed), and DMT(dashed-dotted) 
models. 
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Figure 12: Unpolarized response functions for pir° at 6*£™ = 170° from JLab Hall A. 
The solid line corresponds to a fit to the data using the MAID03 implementation. 
The dashed line corresponds to the SAID solution. 
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Eqs.(35)-(40), the differential cross section in Eq.(17) can then be written as 
,2 2 

= Yl A e p e( cos0 *) + E B e p e( c osd*)]cos(f)M + [C 2 P 2 {cos6*)]cos2(j)M (110) 

The coefficients of the above equation are related to |Mi + | 2 and its projection onto 



the other s- and p-wave multipoles E\ + , S\ + , Mi_, E +, Sq+ ■ 

\M 1+ \ 2 = A /2, (111) 

Re(E 1+ M* + ) = (A 2 - 2C 2 /3)/8, (112) 

Re(M!-M; + ) = -[A 2 + 2(A + C 2 /3)}/8, (113) 

Re(E 0+ M* + ) = Ai/2, (114) 

Re(S 0+ M* + ) = B u (115) 

Re{S 1+ M* + ) = B 2 /6. (116) 



The partial wave coefficents of Eq.(llO) are determined by fitting the differential 
cross sections data such as those displayed in Figs. 9-11. From the relations 
Eqs.(lll)-(116), one then obtains the M 1+ , E 1+ and S\ + amplitudes for determin- 
ing the 7 TV — > A form factors through Eqs.(99)-(101). 

The results from using the above procedure must be corrected for the systematic 
errors due to the truncation of higher multipoles. This can be accomplished by 
calculating the effects of higher partial waves using a realistic paramctrization of 
the higher mass resonances and a realistic model for the background amplitudes. 
For not too large Q 2 values, this method results in reliable multipoles. 

With the above largely model-independent procedure, results for G* M up to 
Q 2 = 6 GeV 2 and the ratios Rem and Rsm, defined in Eqs.(102) and (103), up to 
Q 2 = 4 GeV 2 have been obtained at JLab and are compared with various theoretical 
predictions in Figs. 13 and 14. We now explain how the displayed results from 
SL, MAID, and DMT models are obtained. Within the MAID model, the Q 2 - 
dependence of the 7 V — > N* transition strengths A a of Walker's parameterization 
Eq.(64) is determined from fitting the differential cross section data. The resulting 
multipolc amplitudes are then used to extract the 7./V — > A form factors by using 
Eqs.(99)-(101). On the other hand, within the SL and DMT dynamical models, 
the parameters of the bare quantities Gm{Q 2 ), Ge(Q 2 ) and Gc(Q 2 ), defined by 
Eq.(93), are adjusted to fit the data. The dressed form factors of the SL model are 
then predicted by using Eq.(84) to calculate the meson cloud effect. As shown in 
Ref. 68 , at the A mass W = 1.232 GeV this procedure is equivalent to that based 
on Eqs.(99)-(101). The parameters of these three models have been determined by 
using the data up to Q = 4 GeV 2 . The results at Q 2 > 4 GeV 2 are their predictions. 

In Fig. 13, we see that the theoretical predictions of G* M (Q 2 ) at Q 2 > 4 (GeV/c) 2 
from SL , MAID, and DMT models agree well with the new data from Jlab. The 
prediction by Stoler 119 , which is based on a PQCD-motivated model, is also dis- 
played there for comparison. It also agrees well with the data at relatively high 
Q 2 . 
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The dotted curve in Fig. 13 is obtained from setting the pion cloud effect, defined 
by Eq.(84) and illutrated in Fig. 7, to zero within the Sato-Lee model. We see that 
the pion cloud effect is very large at low Q 2 , but becomes much smaller at high Q 2 . 
Clearly this Q 2 dependence plays an important role in getting the agreement with 
the data up to Q 2 —6 GeV 2 . It will be interesting to see whether the predicted pion 
cloud effect will agree with the data at even higher Q 2 . 




Figure 13: The data of magnetic form factor G* M for the 7-/V — > A(1232) transition are 
compared to various models. Results from old single arm electron scattering experiments 
are labeled "inclusive" . All other results have been obtained from a multipole expansion 
of exclusive ir° production from protons. 

In the upper part of Fig. 14, we see that the preliminary CLAS data for Rem 
at low Q 2 < 0.2 (GeV) 2 139 are in good agreement with the predictions from the 
SL and DMT models. On the other hand, the new Jlab data for the ratio Rsm 
(lower part of Fig. 14) in the low Q 2 < 0.2 GeV 2 region prefers the prediction from 
the SL model. The data points at Q 2 = 0.127 GeV 2 from MAMI and Bates have 
a larger magnitude for Rsm- These data points were used by DMT in fixing their 
parameterization for Gc(Q 2 )- It should be noted that the points from MAMI and 
Bates are not the result of an independent multipole fit, but are from data sets 
with more limited angle coverage fitted to the MAID parametrization. One of the 
data sets from MIT-Bates is shown in Fig. 15. Clearly, only very limited angles are 
covered. Nevertheless, these data are very useful in revealing the non-zero values 
of the Ge and Gc form factors within the dynamical model, as illustrated in the 
difference between the solid and dotted curves. 

We now note that the ratios Rem and Rsm calculated from the dynamical 
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Figure 14: Ratios Rem (denoted as E\ + /Mi + ) and Rsm (denoted as S1+/M1+) for the 
7iV -> A(1232) transition. These two ratios are related to the E\ 12 , Sf( 2 , and Mf{ 2 
multipole amplituds of 7*iV — » tvN, as denned in Eqs.(102)-(103). Preliminary data from 
CLAS at low Q 2 are also included. 
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Figure 15: The data of Alt from MIT-Bates are compared with the results from 
Sato-Lee model. The dashed curve is obtained from setting the Gc and Ge of the 
77V — * A form factor to zero. 
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Figure 16: The imaginary parts (Im) of the E^J 2 and Sf+ 2 calculated from the 
Sato-Lee model. The dotted curves are obtained from setting the pion cloud effect 
to zero. 

models are very much related to the predicted pion cloud effects. These are illus- 
trated in Fig. 16 from the SL model. We see that the pion cloud effect can strongly 
enhence the the E^ 2 and SJ^ 2 amplitudes of 7* TV — > 7riV at low Q 2 . As defined 
in Eqs.(99)-(101), these two amplitudes are related to the Ge(Q 2 ) and Gc(Q 2 ) of 
77V — > A transition. The non-trivial pion cloud effects shown in Fig. 16 are clearly 
verified by the JLab data, as seen in Fig. 14. 

To further improve the determination of the 'fN — > A form factors, data of 
polarization observables must be included in the theoretical analyses. Measurements 
using polarized electron beam and/or a longitudinally polarized hydrogen target 137 
have yielded data of double spin beam-target asymmetry A et and target asymmetry 
A t . Samples of asymmetry data from CLAS are shown in Fig. 17. The double 
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Figure 17: A et and A t from JLab at W=1.225 GeV, Q 2 = QAQGeV 2 are compared 
with the predictions from SL, MAID, and DMT models. The left panel is integrated 
in a range <f> = (144°, 160°). The right panel is integrated in a range cos#* = 
(-0.6,-0.8). 

polarization asymmetry A et is largely given by the well determined M\+ multipolc 
and is well described by all models. However, significant differences can be seen in 
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the At asymmetry which is sensitive to interferences between the non-resonant and 
resonant amplitudes. The discrepancy in the model descriptions can be attributed 
to their different treatments of the non-resonant amplitudes, as discussed in section 
4. 

Extensive pion electroproduction data in the second and third resonance regions 
have also been obtained using CLAS. Some typical results are shown in Figs. 18. 
Here we see that the displayed theoretical predictions do not agree well with the 
data at W > 1.4 GeV. This is not surprising since the parameters of these single- 
channel models arc fixed by mainly fitting the data at W < 1.4 GeV. Recently, 
fits to these higher W data have been achieved by using a single-channel K-matrix 
model and fixed-t dispersion relations. In these analyses pn and mr + data are 
fitted simultaneously using unpolarized cross section data as well as beam spin 
response function results. It has been found that these very different approaches 
give consistent results, e.g. in the analyses of Aznauryan et al. 53 . It indicates 
that the model-dependence may be relatively small. Nevertheless, the extracted 
resonance parameters must be taken with caution before a rigorous investigation of 
the coupled-channel effect has been carried out. Progress in this direction is being 
made 76 . The results of these fits are discussed in section 5.3. 
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Figure 18: The data of a LT i for y*p — > pty° in the 2 nd resonance region are compared 
with MAID (dashed) and the DMT (dashed-dotted) predictions. 
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5.2. Photoproduction and electroproduction of rj mesons. 

In contrast to the pion with isospin 1=1, the eta is an isoscalar meson with 
no charged partners. As such it can only couple with nucleons to form 1 = 1/2 
resonances. This makes the production of r/'s from nucleon targets an ideal tool to 
separate isospin | N* resonances from isospin § A* resonances. The total photopro- 
duction cross section, shown in Fig. 19, exhibits a rapid rise just above threshold, 
indicative of a strong s-wave contribution near threshold. This behavior is known 
to be due to the first negative parity nucleon resonance, the Sn (1535)iVi/2- , which 
couples with approximately 55% to the Nij channel 141 . The nearby £> 13 (1520) has 
a branching ratio of much less than 1% to this channel? The next higher mass nu- 
cleon resonance with a significant Nrj coupling is the Pn(1710), nearly 200 MeV/c 2 
higher in mass. This fact makes the production of 77's from nucleon targets the reac- 
tion of choice for detailed studies of the electromagnetic transition from the ground 
state to the 5n(1535). The Nr) channel effectively isolates this state from other 
nearby resonances, similar to the A(1232) which is well separated from higher mass 
resonances in the pir° channel. In distinction to the A(1232), whose electromag- 
netic transition form fcators drop rapidly with increasing photon virtuality Q 2 , the 
Sii(1535) remains a prominent resonance even at the highest Q 2 that are currently 
accessible. 

In the following subsections we discuss the status of the electromagnetic produc- 
tion of 77's from nucleons, and analyses to extract the photocoupling hclicity ampli- 
tudes for the 77V — * 5n(1535) transition, and their Q 2 evolution. We finally com- 
pare the results with model predictions. Table 4 gives an overview of the kinematics 
covered in recent 77 production measurements I42,i43,i45,i46,i47,i48,i5i,i52,i53,i54_ 

Tabic 4. Summary of rj production data 



Reaction 



Observable 



W range 
(GeV) 



Q range 
(GeV 2 ) 



Lab. 



7(n/p) -> (n/p)r) 



IP ■ 
IP ' 
IP- 
ep - 



•pr\ 
■ PV 
■PV 
eprj 



da/dSl 

da/dU 

da/d^l 

dcr/dQ 

dtj n /da p 

£ 

T 

E 

<?LT, 

<?LT, &TT, 



ax - 



ea L 

tCF L 



< 2.0 

< 1.7 

< 2.3 

< 1.53 

< 2.3 

< 2.3 

< 2.3 

< 1.53 

< 2.2 

< 2.2 



2.8 ■ 
0.3 



4.0 
4.0 



JLab-CLAS 

GRAAL 

ELSA-CB 

MAMI-TAPS 

GRAAL 

GRAAL 

ELSA 

MAMI-A2 

JLab-Hall-C 

JLab-CLAS 



"Even though the Di3(1520) coupling to Nr) is very small, its close proximity to the 511(1535) 
causes large interferences with the dominant Eo+ transition amplitude of the S11. This in turn 
allows a precise determination of the Di3(1520) — » Nr] branching ratio. 
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E,(MeV) 

Figure 19: Total cross section for n photoproduction from protons. Data from MAMI 
(open circles), GRAAL (full circles), and CLAS (open squares) are shown. The curves are 
the fits explained in section 5.2.4. 



5.2.1. r] photoproduction from protons 

With the new measurements in recent years, the data base for rj photoproduction 
reaction has been improved tremendously. The differential cross section data now 
cover the prj mass range up to W = 2.3 GeV, and are available for most of the 
angular range in the hadronic center-of-mass system. Some of these data are shown 
Figure 20. The data from all three experiments agree well. In the mass region of the 
5n(1535) resonance the angular distributions are nearly flat, indicating dominant 
s-wave components with only slight indications of higher partial wave contributions. 
In the mass region above 1.750 GeV, the angular distributions become increasingly 
forward-peaked, indicating significant non-resonant behavior presumably due to t- 
channel processes. 

GRAAL has measured the beam asymmetries using laser light backscattered 
from the 6 GeV electrons to generate high-energy linearly polarized photons. The 
beam asymmetry E 7 is defined as 

s 1 f^ = o°)-5n(^ = 90°) (U7) 

1 p 7 ^(0 = 0°) + ^(0 = 90°) 

where P 7 is the photon polarization, <f> is the azimuthal angle between the plane 
defined by the linear photon polarization and the hadronic plane defined by the 
photon beam and the prj final state. The measured beam asymmetries are shown 
in Fig. 21. Just above n threshold and at the resonance position E 7 shows a 
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Figure 20: Differential cross section for -yp — > p?y from TAPS, GRAAL, and CLAS are 
shown for photon energies near threshold to 1.875 GeV. The curves are the fits explained 
in section 5.2.4. The solid (dotted) line represents the UIM (DR) analysis. 



symmetric angular distribution, approximately following a sin 2 9* behavior, while 
at higher energies the asymmetry is more forward peaked. The sin 2 9* behavior 
near the resonance pole at the lower energies is prominent in the data, and is also 
reflected in the model descriptions included in Fig. 21. 

Asymmetries have also been measured with a transversely polarized proton tar- 
get at ELSA. The target asymmetry T is given by: 

rp _ 1 gnCt) ~ rffl(i) (118) 

*rft(t) + &(D" 

The arrows indicate the direction of the proton polarization relative to the hadronic 
plane. Results are shown in the right panel in Fig. 21. 

5.2.2. "fp — > Sii(1535) parameters extracted from gobal fit to r] photoproduction 

The Sii(1535) has long been known as a strong nucleon resonance with a large 
branching ratio to the Nrj and Nit channels. However, there have been assertions 
that the strong enhancement near this mass is not due to the excitation of a reso- 
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Figure 21: The fits to the 7p — > pr\ data using the Unitary Isobar Model (solid curves) 
and Dispersion Relation (short dashed curves)approaches. Left panel: Beam asymmetry 
measured at GRAAL. The dotted curve in the 1050 MeV panel shows the fit without the 
fi5(1680) resonance. Right panel: Target asymmetry measured at ELSA. 



nance 155 . One model implies that the state is dominantly a dynamically generated 
YjK resonance 156 . On the other hand, recent Lattice QCD calculations 157 ' 158 show 
that there is a strong 3-quark state at this mass with the spin-parity J p = \ , in- 
dicating that the Sii(1535) is indeed an excited state of the nucleon. 

The photocoupling amplitudes and their Q 2 dependence are powerful tools in 
determining the internal resonance structure and will help solve this controversy. 
For the purpose of this article we consider the Sii(1535) as a baryon resonance with 
well defined quantum numbers. In the following we discuss results obtained in a 
global analysis of all observables in r\ photoproduction with the goal to extract the 
photocouplings amplitudes for resonances coupling to the pr\ channel, especially the 
511(1535). 

A number of analyses have been performed on the pi] photoproduction chan- 
nel 53 > 159 > 161 . Here, we describe the recent global analysis by Aznauryan 53 as it 
allows to also assess the model-dependence of the results. Differential cross sec- 
tions 152 ■ 142 were included as well as polarized beam asymmetries 153 , and polarized 
target asymmetries 149 . All established N* resonances above the Ni] threshold were 
included, i.e. Sn(1535), L>i 3 (1520), Su(1650), L>i 5 (1675), Fi 5 (1680), L>i 3 (1700), 
Pn(1710), Pi3(1720). The cross section data are fitted for photon energies up to 
2 GeV, corresponding to invariant masses in the range W = 1.49 - 2.15 GeV, i.e. 
covering the entire resonance region. The polarization data cover only the range up 
to W — 1.7 GeV. Figure 20 and Fig. 21 show samples of the fit to the cross sections 
and asymmetry data. The Unitary Isobar Model (UIM) and Dispersion-relations 
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(DR) approaches give consistent results for the Sn(1535), Z?i3(1520), and Fi 5 (1680) 
resonances. The first result is the confirmation of a large photocoupling amplitude 
for the 511(1535) which is determined with good precision. The results are sum- 
marized in table 5. They are compared with the range given by the PDG 24 . The 



Table 5. Sii(1535) photocoupling from gobal fit in units (10 3 GeV 



Resonance 


Mass (MeV) 


T (MeV) 


A l/2 


Model 


5n(1535) 


1527 


142 


96 


Isobar model 




1542 


195 


119 


Dispersion relations 




1520-1555 


100-200 


60 - 120 


PDG estimate 



new results are within the upper part of the range given by the PDG. The lower 
range in the PDG value comes from an analysis of pion production data by the 
George Washington University (GWU) group 150 . We note here that the results 
from the global fits are also in good agreement with a combined analysis of 7r and 
r\ electroproduction data. We will discuss this in section 5.3. 

From the fit to the differential cross sections one can then also extract the total 
photo absorption cross section for -q production. The fit results are compared with 
the experimental data in Fig. 19. All three experiments agree well in the region 
where the 5*11(1535) resonance dominates, while there is a discrepancy near 1100 
MeV photon energy. Since the angular distributions agree well, this discrepancy 
must be entirely due to different models used for the extrapolation into the un- 
measured angular regions. This emphasizes the importance of measureing complete 
angular distributions which are now available 143 . 

The global analysis incorporates also the beam asymmetry in the fit. To illus- 
trate the sensitivity of X to contributions from the Di3(1520) we express S 7 in the 
approximation that only S- waves, P-waves, and D- waves with spin J < § contribute 



as 160 : 



_ 3sin 2 gife[g * + (£ 2 -+M 2 -)] 

T ~ c 12 \ LLy ) 



Table 6. Summary of r\ photoproduction gobal fit results. The uncertainty in /3^jv reflects the 
model dependence in using the unitary isobar model and dispersion relation approach. 



Resonance 


Mass (MeV) 


T(MeV) 


P v n(%) 




£>l 3 (1520) 


1520 


120 


0.05 ± 0.02 


50 - 60 


Fi 6 (1680) 


1675 


130 


0.15 ± 0.03 


60 - 70 



This expression can be fitted to the measured beam asymmetry S 7 . Using Eq + 
from fits to the cross section data, the multipoles E%_, M%_ for the Di3(1520) can 
then be determined. Since the corresponding pion multipoles E%_, M£_ are known 
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with high precision from pion production, the branching ratio f3 n N can be extracted. 
The analysis also allows to extract the Nrj branching ratio for the i*i 5 (i680) by an- 
alyzing the forward-backward asymmetry in S 7 seen in Fig. 21 at E 1 w 1 GeV. The 
dotted curve in the figure for 1050 MeV (left panel) shows the fit when the small 
Fi5(1680) amplitudes are turned off. Clearly, the ( interference effects strongly 
enhance this contribution. The results for the Di3(1520) and Fi 5 (1680) are sum- 
marized in table 6. Both results represent significantly improved values for the 
branching ratios. 
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Figure 22: Samples of differential cross section for j*p — ► p»7 from CLAS at fixed Q 2 = 
0.8 GeV 2 and different W and <j) values. The shaded bands indicate the fit results using 
9 Legendre polynomials, as explained in section 5.2.5. The dashed curve represents the 
predictions of jyMAID. 



5.2.3. Eta electroproduction 

Eta electroproduction experiments have focussed on the Q 2 evolution of the 
«Sii(1535) transverse photocouplings amplitude Ai/ 2 {Q 2 )- Experiments at DESY 162 ' 
and Bonn 164 > 165 found a very slow falloff with Q 2 . Recent experiments at Jefferson 
Lab 14 M5i,i54 have studied this behavior in detail with high statistics, and also ex- 
tended the kinematics range. Figure 22 shows samples of differential cross sections 
measured with CLAS 154 . Even at the peak of the «Sii(1535) resonance the angular 
distributions are not completely flat indicating that higher partial waves are present 
in addition to the dominant S-wave. Figure 23 shows samples of total cross sections 
at fixed Q 2 . In contrast to the A(1232) which rapidly drops with Q 2 , the 5n(1535) 
remains prominent even at the highest Q 2 . 

Most of the published results on the Q 2 dependence of the 5n(1535) transition 
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Figure 23: Left panel: Total cross section for y*p — > pr\ vs W, at various photon virtualities 
Q 2 . The right panel shows partial wave fit parameters vs W compared to predictions of 
77MAID. All data are from CLAS. 



amplitude have been obtained in single resonance fits. This has been justified with 
the dominant contributions of the 511(1535) to the 7*p — > prj cross section. It 
is, however, not a fully satisfactory solution, as higher mass states that couple to 
Nr] may also contribute in the lower mass region. The results have to be taken 
with caution. The differential cross sections are fitted to the expression Eq.(17). 
The dependence on the 77 scattering angle (9*) can be examined by expanding each 
component of the differential cross section in terms of Legendre polynomials: 

00 

J2^Pe(cos6* v ) (120) 

e=o 
00 

J2 B * p t (cos e;) (121) 

£=1 
00 

J2CzP? (cos 9* n ) (122) 

If the expansion is limited to I = 2, only the coefficients Aq, A\, A2, B\, B2 and 
C2 are retained. Results from the fit at fixed Q 2 are shown in the right panel of Fig. 
23. Strong variations of A\ and B\ are seen in the W range from 1.6 to 1.7 GeV, 
indicating large interference effects involving s- and p-waves. Possible p-wave can- 
didates are the Pn(1710) and Pi3(1720) states. A is mostly due to the 5n(1535) 
resonance, and is the by far largest amplitude. The longitudinal and transverse re- 
sponse functions cannot be separated in this analysis. In earlier experiments 165 > 163 
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the longitudinal and transverse cross sections were separately determined at some 
fixed Q 2 values, showing that Ao is dominated by the transverse amplitude A 2 , 2 
in the Q 2 range of this study. The combined analysis of 7r and rj electroproduction 
data, to be discussed in the next section, also finds small longitudinal contribution 
to Aq. Assuming — 0, |j4i/2| can be computed from Aq. Figure 24 shows a 
compilation of results for the 7p — ► 5n(1535) photocouplings hclicity amplitude 
Ax/2{Q 2 )- The slow falloff with Q 2 confirms the unusually hard transition form 
factor that persists to the highest measured values of Q 2 . The solid and dotted 
curves are the prediction of Close and Li 177 , and of Giannini, Santopinto and Vas- 
sallo 175 . It should be noted that the absolute normalization of the data displayed 




Figure 24: Q 2 evolution of the -yp — ► 5n(1535) photocouplings amplitude A 1 / 2 - The 
full symbols are the most recent CLAS results. The open diamonds are the JLab data 
from Hall C, and the open circles are from previous experiments. The solid curve is from 
Ref.(177), and the dotted curve is from Ref.(175). 

in Fig. 24 is uncertain to the extent that the branching ratio (3nt)(Sh) — 0.55 and a 
total widths of 150 MeV have been used in extracting Ai/ 2 - The Review of Particle 
Properties 2002 allows a large range of 0.30 - 0.55 for the branching ratio. However, 
the recent analysis of Armstrong et al. 141 , gives a value of /3at 7? (S'ii) ~ 0.55. The 
use of this value is consistent with the values Pnti(Sh) = 0.55 and (3ntt(Sh) = 0.4 
used in the combined analysis of it and rj electroproduction which is the subject of 
next section. 

5.3. Combined analysis of n and rj electroproduction data 

The large amount of data taken by the CLAS detector allows simultaneous 
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measurements of cross sections and polarization observables for several channels, e.g. 
pir°, mr + , prj. Also, the large acceptance provides complete angular distributions, 
including the full azimuthal dependence. Use of a highly polarized electron beam 
provides data on the helicity-dependent response function <jlt' covering the full 
angle range. Results of a MAID and DMT analysis of the channel p-K° have recently 
been reported 138 . Combined analyses of these data in a multi-channel global fit 
provides much more stringent constraints on resonance parameters than single- 
channel analyses can. The full set of data taken with a hydrogen target have been 
analyzed within the unitary isobar model 52 and the dispersion relation approach 53 
described in section 4. The data on <tlt' are especially sensitive to small resonance 
contributions in a large non-resonant background. The sensitivity is the result of 
the interference term that mixes real and imaginary amplitudes 

a lt' ~ Im(L) ■ Re(T) + Im(T) ■ Re(L) , (123) 

where L and T represent the longitudinal and transverse amplitudes, respectively. 
Figure 25 shows the the sensivity of ctlt' to the Pn(1440) multipoles Mi_ and 
Both channels show sensitivity to changes in the multipoles, however, the effect is 
much larger in the nir + channel. This is due to a combination of two factors, the 
stronger coupling of I = 1/2 nuclcon states to the mr + channel, and the larger 
background terms contributing to this channel. Using pion and eta production 
to study the same resonances, e.g. the 5n(1535) allows tests of possible flavor- 
dependence of the results which could be significant if rescattering effects, such as 
present in the A(1232) region, play an important role also for higher mass states. 
The use of two conceptually very different approaches also allow to estimate the 
model-dependence of the resulting amplitudes. The results for the Pn(1440) are 
shown in Fig. 26. For the first time a consistent trend is emerging: The magnitude 
of A 1 / 2 {Q 2 ) drops rapidly for Q 2 > 0, with a sign change near Q 2 = 0.5 GeV 2 . 
There is also a strong longitudinal coupling. Bold, solid, dashed and dot-dashed 
lines are from various calculations 171,173 ' 174,176 . Nonrelativistic quark models 172 ' 177 
predict large negative A 1 / 2 (Q 2 ) in the entire Q 2 range, and do not describe the 
data. The hybrid model 171 describes the fast drop of | / 2 j qualitatively, but has 
no sign change, and predicts Si/2(Q 2 ) = 0, while the data show a sizeable Si/ 2 
amplitude. The relativistic models of Capstick and Keister 173 and of Cardarelli 
and Simula 174 predict the sign change for A 1 / 2 (Q 2 ) but show a much faster rise 
than is observed. The magnitudes and trends of both amplitudes are well described 
by a model that describes the Pn(1440) with a small quark core and a qq cloud 176 . 
In this model, the low Q 2 behavior is entirely due to the qq contribution while the 
qqq core defines the high Q 2 behavior. The results for the Sn(1535) are depicted 
in Fig. 27. They show consistent results for A\i 2 in the prj and the Nit channels for 
the Unitary Isobar Model (UIM) and Dispersion-Relation (DR) analyses. Also, for 
the first time, stable results for the longitudinal coupling Si/ 2 have been obtained. 
The solid and dotted curves in Fig. 27 represent quark model calculations using a 
harmonic oscillator potential 177 and a hypercentral potential 175 , respectively. 
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Figure 25: Sensitivity of CLAS data of a hT , to changes in multipoles (dashed) and 

(dotted) for mv + channel (left) and pn° channel (right). Solid line shows best fit 
using the Unitary Isobar Model. The -0.5 in the legend refers to shifting the Breit-Wigner 
amplitude by -0.5 /ifc 1 / 2 . 
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Figure 26: Results for the Pn(1440) amplitudes A p 1/2 and S^ 2 . CLAS points include 
model error. The Hall A point shows MAID03 fit model error. The curves are described 
in the text. 
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Figure 27: Q 2 evolution of the Sii(1535) photocoupling amplitudes Ai/2(Q 2 ) (left), and 
S\/2{Q 2 ) (right). Cross section data from pr/ y pit , and nir + have been used, as well 

as polarized beam response function cr£y/ and <rJ T / . The shaded band indicates the 
uncertainties seen in previous analysis using prj cross sections data. Solid and dotted lines 
are from quark model calculations described in the text. 
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Figure 28: Q evolution of the Di3(1520) photocoupling amplitudes A1/2, A3/2, and S1/2. 
Cross section data from pn° , and n-K + have been used, as well as polarized beam response 

+ o 

function <t£ t , and (t£ t /. The shaded bands indicate the uncertainties seen in previous 
analysis using mostly pn° cross section data. The theoretical curves are explained in the 
text. 
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The results of the global fit for Ai/ 2 , A^/ 2 and S\/ 2 for the Di3(1520) are shown 
in Fig. 28. Both the UIM and DR analyses give consistent results. To summarize, 
the inclusion of polarization observables in addition to the differential cross section 
into a global analysis results in a less model-dependent description of 7r and rj 
photoproduction and electroproduction processes in the resonance region. There 
are several noteworthy results: First, consistent results are obtained for the mass 
of the «Sii(1535) for the 77 and it channel, both for photo- and electroproduction. 
The mass is in the range M(Sn) = 1531 ± 5 MeV. Second, the discrepancy between 
r\ and 7r photoproduction results for A X j 2 amplitude seems to have been resolved. 
The analysis of electroproduction data gives also good agreement for the pr] and 
Nit channels and in both UIM and DR approaches. Third, the Q 2 evolution of 
the A 1 / 2 and Si/ 2 amplitudes for the Pn(1440) are consistent with the predictions 
of a meson cloud model. This is in line with what has been found earlier for 
the 7VA(1232) transition, that meson cloud effects can be sizeable for some of the 
resonance transitions. 

5.3.1. Analysis of Resonance Transitions in the Single Quark Transition Model 

Properties of nucleon resonances such as mass, spin-parity, and flavor fit well 
into the representation of the 5/7(6) £g> 0(3) symmetry group, which describes the 
spin-flavor and orbital wave functions of the 3-quark system. This symmetry group 
leads to supermultiplcts of baryon states with the same orbital angular momentum 
L of the 3-quark system, and degenerate energy levels. Within a supermultiplct 
the quark spins are aligned to form a total quark spin s, with s = |, §, which 
combines with the orbital angular momentum L to form the total angular momen- 
tum J = L + s. A large number of explicit dynamical quark models have been 
developed to describe the electromagnetic transitions between the nucleon ground 
state and its excited states 169 . 172 ,i73,i75^ Measurement of resonance transitions 
and the dependence on the distance scales, given by the virtuality of Q 2 of the pho- 
ton, provides information on the nucleon wave function. In order to compute the 
transition, assumptions on the 3-quark potential and the quark-quark interactions 
have to be made. These are then tested by predicting photocoupling hclicity ampli- 
tudes which can then be confronted with experimental data. Algebraic relations 
have been derived for resonance transitions assuming the transition only affects a 
single quark in the nucleon. The parameters in these algebraic equations can be 
determined from the experimental analysis 178 . Based on the symmetry properties 
of the Single Quark Transition Model (SQTM), predictions for a large number of 
resonances belonging to the same SU(6) (g> (9(3)supermultiplet can be made. The 
fundamentals of the SQTM are described in references 179 > 180 ; where the symmetry 
properties have been discussed for the transitions from the ground state nucleon 
[56, + ] to the [70, 1~] and the [56,2+] supermultiplets. The [70, 1~] contains states 
which arc prominent in electromagnetic excitations, and it is the only supermulti- 
plet for which sufficient data on resonance couplings of two states are available to 
extract the SQTM amplitudes and test predictions for other states. The coupling 
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Q 2 (GeV 2 ) 

Figure 29: Single quark transition amplitudes A, B, C as functions of Q 2 . The 
amplitudes have been extracted from the recent JLab, MAMI, and GRAAL data 
on the 511(1535) and world data on £>i3(1520) photocoupling amplitudes. The 
shaded band is a parameterization of the data from Ref. 178. 



Table 7: Helicity amplitudes for the electromagnetic transition from the ground 
state [56,0+] to the [70, 1~] multiplet as a function of the SQTM amplitudes. 9 
is the mixing angle relating two J p = | states with s 3q — | and s 3q = |. 

There is also a small mixing angle for the two | states resulting in the physical 
states £>i3(1520) and Di3(1700). We have not included the latter mixing angle in 
the table. Note that the excitation of the Z?i3(1700) from a proton target is only 
possible because of SU(6) symmetry breaking leading to the mixing with the lower 
mass Di 3 (1520). 



State Proton target Neutron target 



Su(1535) 


A + 


U A + B 


- c) cos e 


A° 

A l/2 - 


-I(A+lB-IC) 


Di 3 (1520) 


A+ 




2B-C) 


A° 

A l/2 - 


~^ A ~ 2B '^ 




A + 

A 3/2 ~ 


-±=(A + 

2V6 y 


C) 


4° 

A 3/2 ~ 


^ 3A -^ 


Sn(1650) 


A + 

A l/2 ~ 


\{A + B 


- C)sin0 


4° 

A l/2 ~ 




Di 3 (1700) 


4 + 

A l/2 - 







A° 

A l/2 - 


I57f( B - 4C ) 




A + 

A 3/2 - 







A° 

^3/2 _ 




Di 5 (1675) 


A + 

A l/2 ~ 







A° 

A l/2 - 






A + 

A 3/2 - 







A° 

A 3/2 — 


-ly/liB + C) 


D 33 (1700) 


A l/2 = 


^ A - 


2B-C) 


same 






^3/2 _ 


C) 


same 




531(1620) 


A l/2 = 


U* A - 


B + C) 


same 
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of the electromagnetic current is considered for the transverse photon component, 
and the quarks in the nucleon are assumed to interact freely with the photon. In 
such a model the quark transverse current can be written in general as a sum of 
four terms 179,180,181. 

J+ = AL+ + Ba+L z + Ca z L+ + Da~L + L + , (124) 

where a is the quark Pauli spin operator, and the terms with A, B, C, D in front 
operate on the quark spatial wave function changing the component of orbital an- 
gular momentum along the direction of the momentum transfer (z- axis). The A 
term corresponds to a quark orbit flip with AL Z = +1, term B to a quark spin 
flip with AL Z — 0, the C and D terms correspond to simultaneous quark orbit and 
quark spin flip with orbital angular momentum flips of AL Z — +1 and AL Z = +2, 
respectively. For the transition from the [56, + ] to the [70, 1~] supermultiplet with 
L = 1, only A, B, and C are allowed. The relationship between the A, B, C am- 
plitudes and the usual helicity photocoupling amplitudes A 1 / 2 and A 3 / 2 is listed 
in table 7. Using the extracted photocoupling amplitudes from the 5 l n(1535) and 
the L>i 3 (1520), the A, B, C amplitudes for the 7 + [56,0+] -> [70, 1"] have been 
extracted 178 . The results are shown in Fig. 29. Knowledge of the 3 amplitudes 
and of two mixing angles for the transition to the [70, 1~] allows predictions for 16 
amplitudes of states belonging to the same supermultiplet. If they can be confirmed 
for some of the amplitudes, one then has a measure of the degree to which electro- 
magnetic transitions of nucleon resonance are dominated by single quark transitions 
at the photon point (Q 2 = 0) and, using electroproduction data, examine if and 
how this is changing as a function of the distance scale at increasing photon virtu- 
ally. The SQTM predictions for the proton and neutron amplitudes are shown in 
Fig. 30. There is remarkable agreement between the predictions and the data at 
the photon point. For electroproduction, there is good agreement where consistent 
data sets are available, i.e. for the Sii(1650). Much improved electroproduction 
data are needed for more definite conclusions. Most of the states belonging to the 
[70, 1~] supermultiplet with masses near 1700 MeV couple strongly to Ntttt chan- 
nels. Studies of these channels require use of large acceptance detectors and new, 
sophisticated analysis techniques. Progress made in these areas will be discussed 
in the next section. There are similar relations for the transition from the nu- 
cleon ground state to the members of the [56, 2 + ] supermultiplet. In this case four 
SQTM amplitudes can contribute. Unfortunately, the only state for which the two 
transverse photocoupling amplitudes have been measured in electroproduction is 
the Fi 5 (1680). This is insufficient to extract the four SQTM amplitudes. 

5.4. Two-pion production. 

Two-pion channels dominante the electromagnetic meson production cross sec- 
tions in the second and third resonance regions where we hope to resolve the missing 
resonance problem 169 and ultimately determine what basic symmetry group 206 is 
underlying the baryon spectrum. Thus, a detailed understanding of two-pion pro- 
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Figure 30: Single quark model transition prediction for the [70, 1~] multiplct. The 
SQTM predictions are shown by the shaded band in comparison with the experi- 
mental data. At Q 2 = the full circle is the Particle Data Group estimate. For 
Q 2 > 0, measurements from JLab, Bonn, DESY, and NINA in rj and it clectropro- 
duction are shown. For the 5 l n(1535), the results of an analysis of the world data 
in 77-electroproduction presented in Rcf. 151 are also included. The superscript 
refers to neutron states. 



•55 



duction is very important in the N* study, and has been pursued very actively in 
recent years. Very extensive two-pion production data have now been accumulated 
at JLab, MAMI and ELSA, but have not been fully anaylzed and understood the- 
oretically. Here we will mainly report on the status of the data and describe some 
very preliminary attempts to identify N* states. 

The study of Nitty channels requires the use of detectors with nearly 4n solid 
angle coverage for charged or neutral particle detection. Several such detectors have 
been in operation for a number of years, and have generated large data sets for the 
following reactions 

7P — > ptt tt , (125) 
7P — ► p7r + 7r~ , (126) 
ep — > epn + ir~ . (127) 

The two-pion channels in the above processes can be projected onto various 
isobar channels which are more useful in identifying the nucleon resonances from the 
data. The pn + -K^ final state is sensitive to the A ++ ir~ isobar channel which could 
have large sensitivity to resonance decays but has also very strong contributions 
from non-resonant mechanisms. The pir°ir° final state has the advantage of high 
sensitivity to resonance contributions with fewer non-resonant contributions. The 
pir + n~ channel is sensitive to the pp° isobar channel, while the p7r°7r° channel does 
not couple to pp°. Table 8 gives an overview of recent 2-pion production data 
obtained at various laboratories. 

Table 8. Summary of •yp — > pnn production data 

Reaction Observable W range Q A range Lab. 
(GcV) (GcV 2 ) 



7 p — > pTV^TT 



7P — * P 71 " 7r 



atot, M p7r+ , M n+7V - 



<?tot, M„ 



da 



&t0t, M p7r O , ¥,0,0 

atot, 2, M pl o, M n o n o 
event-by-event analysis 



< 2.1 0.65 

< 2.7 

< 2.0 

< 1.9 

< 1.55 

< 2.6 



1.3 

1.5 - 4.0 








CLAS 184 

CLAS 182 
CLAS 186 

GRAAL 188 
MAMI 189 ' 190 
CB-ELSA 183 



5.4.1. Analysis of the data withpn + n final state 

The quality of the recent data with ptt + it~ final state is very high. An example 
is shown in Fig. 31 for the differential cross sections in one-dimensional projections. 
The data show evidence for the formation of the A ++ 7r~, A°7r + , and pp° isobar 
channels. In the real photon case (left panel) the pp° contribution dominates the 
higher invariant mass (W) region. This contribution drops significantly for virtual 
photons (right panel). We also see that the tt" angular distribution is much more 
forward-peaked in the case of real photons, and becomes flatter with increasing Q 2 . 
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This behavior can be qualitatively understood as a consequence of vector meson 
photon duality. With increasing Q 2 the probability of the photon to fluctuate into 
a virtual vector meson is reduced leading to reduced diffraction-like scattering. 




Figure 31: One-dimensional projections for the 7p — > pir + n reaction cross sections. Left: 
Q 2 = 0. Right: Q 2 = 0.95 GeV 2 . Both data set are from CLAS. Curves are explained in 
the text. 




W (GeV) W (GeV) 

Figure 32: Total cross section for photoproduction (left) and electroproduction (right) of 
pir + iv~ at Q 2 — 0.65, 0.95, 1.30GeV 2 (from the top). Data sets are from CLAS. Curves 
are explained in the text. 

Although the high quality data of two-pion production are now available, the 
analysis in terms of extraction of resonance parameters has not been fully developed. 
The theoretical understanding of these complex processes, mainly based on the 
tree-diagram isobar approach 90 , is also very limited and preliminary. Recently, two 
distinctly different approaches have been applied to analyze the photoproduction 
and electroproduction data from JLab 186 ' 191 and CB-ELSA 183 . The first approach 
is to adjust the parameters of an isobar model 191 to fit the fully extracted cross 
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WGeV/c 2 WGeV/c 2 



Figure 33: Preliminary results of a partial wave analysis in the "/p — > pn + n showing the 
m = i spin projections for the well-known Fi5(1680) (left), and evidence for the poorly 
known P 33 ( 1600) (right). 



section and polarization asymmetry data 191 < 192 . The second one is to fit directly 
the unbinned data event-by-event 183 . We describe them in the remaining part of 
this subsection. 

The first approach makes use of knowledge from hadronic production. The 
energy-dependence of non-resonant processes is parameterized, and resonance pho- 
tocouplings and hadronic couplings are fixed, if known, e.g. from single pion pro- 
cesses. Resonances in specific partial waves can be introduced to search for undis- 
covered states. Model parameters are usually fitted to the one-dimensional projec- 
tions of the multi-dimensional differential cross section. Such a model can lead to 
a qualitatively good description of the projected data as shown in Fig. 31. The 
method has been used in the analysis of CLAS electroproduction data 18 . In this 
analysis a significant disprepancy was found near W = 1.7 GeV between the data 
and the resonance parametrizations implemented in the fit model. This discrepancy 
was attributed to either inaccurate hadronic couplings for the well known Pi3(1720) 
resonance determined in the analysis of hadronic experiments, or to an additional 
resonance with J p = | + with either I = ~ or J = |. The discrepancy is best 
visible in the total cross sections for electroproduction, shown in Fig. 32. The dot- 
ted line shows the model predictions using resonance parameters from single pion 
electroproduction and from the analysis of ttN — > Nirir data 24 < 58 * 63 . The solid 
line represents the fit when the hadronic coupling of the Pi3(1720) to Air and N p 
are allowed to vary much beyond the ranges established in the analysis of hadronic 
data. Alternatively, a new state was introduced with hadronic couplings extracted 
from the data. Table 9 summarizes results of the analysis using a single Pi 3 with 
modified hadronic couplings, and a new P/3 state with undetermined isospin while 
keeping the PDG Pi3(1720) hadronic couplings unchanged. In either case, the fit 
requires a resonance with hadronic couplings that are significantly different from 
the ones of the Pi 3 (1720) state listed by PDG. 

The total photoproduction cross section in the left panel of Fig 32 shows a W 
dependence that is very different from the electroproduction data in the right panel. 
In particular, the photoproduction has a much higher nonresonant contribution 
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Table 9. PDG parameters for the Pi3(1720) and parameters resulting from fits to pir + n electro- 
production data. 





Mass (McV) 


T(MeV) 


rwr(%) 


T Np /V{%) 


PDG P13 


1725 ± 20 


114 ± 19 


63 ± 12 


19 ± 9 


PDG 


1650 - 1750 


100 - 200 




70 - 85 


New P I3 


1720 ± 20 


88 ± 17 


41 ± 13 


17 ± 10 



largely due to increased non-resonant p° production at the photon point. Both 
data are, however, consistent with a strong resonance near W — 1.72 GeV in the 
P13 partial wave 185 . The drawback of the approach described above is that when 
fitting one-dimensional projections of cross sections, correlations between the data 
sets are lost. 

The second approach 186 is based on a partial wave formalism starting from the 
T matrix at a given photon energy E: 

T fi {E) = <jm+TT-;T f \T\ 1P ;E> 

= ^ < p7T + 7r~; Tf\a >< a\T\jp: E > 

a 

= J2r(r f )V a (E) , (128) 

a 

where a denotes all intermediate states, and r/ characterizes the final state kine- 
matics. The decay amplitude ip a {rf) =< pn + TT~;Tf\a > is calculated using an 
isobar model for specific decay channel, e.g. A ++ 7r~, A~n + , or pp° . The produc- 
tion amplitude V a =< a\T\^fp; E > is then fitted at fixed energy using an unbinned 
maximum likelihood procedure. This method makes use of all information contained 
in the data, and takes into account all correlations between the variables. 

In this analysis a total of 35 partial waves were included in addition to t-channel 
processes with adjustable parameters. Figure 33 shows intensity distributions in 
different isobar channels, for the |(m = |), and |(to = |) partial waves. Clear 
signals of the i*i 5 (1680) and the P33(1600) are seen, the latter being a not fully 
established 3-star resonance. In the final analysis the energy-dependence is fitted 
to a Brcit-Wigncr form to determine masses and widths of resonant states. This 
method is closer to a model-independent approach, and can directly 'discover' new 
resonances in specific partial waves. 

5.4.2. Description of jp pit ir° in resonance analyses. 

The CB-ELSA collaboration has analysed the pir°ir° final state using a more 
model-dependent version of the partial-wave-analysis described above. Here s- 
channel Breit-Wigner distributions are fitted to the data on an event-by-event basis, 
therefore retaining the correlations in the data. However, the fit is constrained by 
the parametrized energy-dependence of the Brcit-Wigncr function. In Fig. 34 the 
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total cross section for jp — * tt°tt is shown as extracted from the integral over 
all partial waves contributions in comparison to previous data from TABS 189 and 
GRAAL 188 . 
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Figure 34: Total cross section of 7p — > pn°n° reaction. The connected line is from the 
preliminary results of a partial wave analysis. The data are from the measurements with 
TAPS and GRAAL. 

5.5. Kaon production. 

Production of kaons from nucleons has long been recognized as a potentially 
very sensitive tool in the search for excited baryon states 103 . Analyses of the KA 
and K £ channels include the isospin selectivity; the KA final state selects isospin 
i, similar to the Nn channel, while KT, couples to both N* and A* resonances. An 
important tool in resonance studies is the measurement of polarization observables. 
The self-analyzing power of the weak decay A — > p7r _ can be utilized to measure 
the A recoil polarization. To make full use of this unique feature large acceptance 
detectors are needed. 

For a long time, the lack of consistent data sets for KA and KY> production in a 
wide kinematics range has hampered the use of kaon production in the study of non- 
strange baryon resonances. Moreover, the interpretation of these data, mainly for 
charged K + A and K + Yj° channels, in terms of N* excitation is complicated by the 
fact that they may be dominated by the nonresonant particle-exchange processes. 
Another drawback in comparison to Nir and Ntttt is the relatively small cross 
section, and the lack of known strong resonances with a dominant coupling to kaon- 
hyperon channels. This fact makes it more difficult to use strangeness production 
as a tool in the study of excited baryons, and specifically in the search for new 
resonances. 



Most of the available theoretical models 79 - 80 ' 81 - 82 ' 83 - 84 ' 193 ' 194 - 195 for kaon pro- 



duction are based on the tree-diagrams approach, as described in section 4.2. The 
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Figure 35: Left: Angular distributions of K + A photoproduction at fixed photon energy 
Ej. Right: K + T,° angular distributions. Both data sets are from SAPH1R. The lines 
represent Legendre fits to the data. 



validity of these tree-diagram models is questionable, as discussed, for example, in a 
coupled-channel study 74 of Kaon photoproduction. We therefore will mainly focus 
on the status of the data, not on the results from these theoretical models. 

5.5.1. Photoproduction of K +A and K + H 

High statistics data of kaon photoproduction covering the resonance region are 
now available from the SAPHIR 196 and the CLAS 197 collaborations. These new 
data consist of high statistics angular distributions as well as A polarization asym- 
metries, as summarized in table 10. 

There are significant discrepancies between the CLAS data and the published 
SAPHIR 187 data, while the new data from SAPHIR 196 are in much better agree- 
ment with the CLAS data. We therefore disregard the earlier published data from 
SAPHIR. Unfortunately, most of the model calculations have been fitted to the 
published results, and therefore can not be reliably compared to the new data. 

The angular distributions for K + A and K + T, production are shown in Fig. 35. 
We see that the K + A data (left panel) show a strong forward peaking for pho- 
ton energies greater than 1 GeV, indicating the large t-channcl contributions. For 
the K + Yi° channel (right panel) the angular distribution are more symmetric or 
"resonance-like" at low energies , but become soemwhat more forward-peaked at 
energies above 1.3 GeV. 

The high statistics of these data allows, for the first time, to identify the struc- 
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Table 10. Summary of hyperon photo- and electroproduction data 



Observable W range cos 0^ range Experiment 
(GcV) (GeV 2 ) 



7P — » AX da/dQ 
da/dQ 
da/dQ 

■yp -> EE" dcr/dQ 
da/dQ 

jp^K+A,S P 



< 2.15 

< 2.6 

< 2.3 

< 2.15 

< 2.6 

< 2.3 

< 2.3 

< 2.6 

< 2.5 

< 2.15 



-0.95 - +0.95 
-0.95 - +0.95 
-0.85 - +0.85 
-0.95 - +0.95 
-0.95 - +0.95 
-0.85 - +0.85 
-0.85 - +0.85 
-0.95 - +0.95 
-1.0 - +1.0 
-1.0 - +1.0 



SAPHIR 187 
SAPHIR 196 
CLAS 197 



SAPHIR 187 
SAPHIR 196 
CLAS 197 




CLAS 197 
SAPHIR 196 
< 3 CLAS 199 
0.3 - 1.5 CLAS 200 



tures in the differential cross section that hint the interference between the reso- 
nances and the nonresonant backgound. The presence of s-channel resonances is 
particularly evident in the W-dependence of the differential cross section shown in 
Figure 36. At the most forward angles (upper panel), two resonance- like structures 
are visible at W w 1.7GeV, and at W w 1.95 GeV. The structure at 1.7 GeV 
could be accomodated by the known states 5n(1650), Pn(1710), and Pi3(1720), if 
the KA coupling of these states is allowed to vary. From hadronic processes these 
couplings are very poorly known 24 . At intermediate angles (middle panel) the data 
indicate a smoother falloff with W, while at backwards angles (lower panel) another 
resonance-like structure near W w 1.875 GeV emerges, overlapping with the struc- 
ture at the higher mass. These distributions reveal complex processes, indicating 
contributions from more than a single resonance near W = 1.9 GeV. 

Samples of the A polarization measured with CLAS are shown in Figure 37. 
The data show a strong W dependence especially at backward angles (upper panel). 
Comparing the displayed data at two angles, we can conclude that the angle distri- 
butions of the A polarization change sign from largely negative at forward angles 
to positive at backward angles. 

The three theoretical results displayed in Figs. 36 and 37 only describe very 
qualitatively the main features of the data in the region of the nucleon resonance 
region. 

5.5.2. Electroproduction of K + A and K + H. 

Kaon electroproduction is another tool in the study of non-strange nucleon res- 
onances. While the K + A and K + H° photoproduction cross section exhibits a com- 
plex structure of resonant and nonresonant contributions that is difficult to disen- 
tangle, some of the resonance contributions in electroproduction may be enhanced 
at higher Q 2 due to their slower form factor falloff compared to other resonances, 
and compared to the background amplitudes. A significant amount of data has be- 
come available recently from CLAS 198 < 200 . In these experiments the electron beam 
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Figure 36: W-dependence of K + A production for several bins in cosO^. Points with full 
circles are from CLAS. The triangles are older data from SAPHIR. The theoretical curves 
are from Ref.201 (Guidal/Laget/vdH), Ref.194 (KAON-MAID), and Ref.195 (Janssen). 



is polarized, and hence the virtual photon also has a net circular polarization. 

Figure 38 shows samples of the K + A production cross sections integrated over 
either the forward hemisphere (left panel) and backward hemisphere (right panel) 
at fixed Q 2 . The results reveal resonant behavior near W=1.7 GeV and 1.87 GeV 
at large angles while at the forward angles the resonant structures are masked by 
the large non-resonant contributions. The enhancements in the cross section appear 
in the same invariant mass W range as in photoproduction, and are likely due to 
the same resonances contributions. 

The data of A recoil polarization have been obtained in measurements with 
polarized electron beams. The measured total A recoil polarization can be written 
as 

P A = P° ± PeP' (129) 

where P e is the electron beam polarization, P° is the induced polarization which is 
present without beam polarization, and P' is the transferred polarization. Figure 
39 displays the data of the transferred A polarization integrated over all Q 2 for 
three bins in W. The considered P' x , and P' z , are the projections of the polarization 
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Figure 37: A polarization in K + A photoproduction measured with CLAS. The theoretical 
curves are the same as those displayed in Fig. 36. 



vector P' onto the x' and z' axes which are also defined in Fig. 39. The data show 
that the z'-polarization is large and is rising with cos 6^, indicating a t-channcl 
mechanism. On the other hand, the x' polarization is large and remains negative 
throughout the angular range. None of the displayed theoretical results from tree- 
diagram models 81 > 194 > 195 and a Regge model 201 can give an adequate description 
of the data. 

To summarize, production of K + A and K + T.° from protons exhibit evidence 
of s-channel nucleon resonance contributions in the mass range where no A^* or 
A* resonances have been well established. However, resonances are masked by 
large t-channel processes. In order to extract reliable information on contributing 
resonances a better understanding of nonresonant processes is needed. Currently, 
the most important task is to continue experimentally to establish a broad and solid 
base of consistent data in the strangeness sector, including extensive differential 
cross sections, beam and target polarization asymmetries, and polarization transfer 
measurements. A "complete" measurement of all observables which is needed to 
unambiguously extract all helicity amplitudes can be achieved 203 > 204 . This requires 
use of a polarized photon beam and of a polarized target and the measurement of 
the hyperon recoil polarization. Experimental effort in this direction will continue 
with a series of new measurements planned at JLab 144 . On the theoretical side, 
a dynamical coupled-channel approach, such as that described in section 4.6, must 
be developed to interpret the extracted N* parameters. 
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Figure 38: Cross section from CLAS for K + A electroproduction at Q 2 — 0.7GeV 2 inte- 
grated over to the forward hemisphere (left panel) and backward hemisphere (right panel) 
in the center-of-mass angle 9* K . 




Figure 39: Left: Transferred A polarization in electroproduction from CLAS. The data 
have been integrated over all <f> angles. The curves are discussed in the text. Right: 
Coordinate system defining the A polarization projections. The curves are predictions of 
models in Ref. 81 (dotted), Ref. 194 (solid), and Ref. 195 (dot-dashed). 
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5.5.3. Photoproduction and electroproduction of vector mesons. 

The early investigations of photoproduction and electroproduction of vector 
mesons were mainly in the high energy region where the data can be explained 
largely by the diffractive Pomeron-exchange mechanism. For the study of nucleon 
resonances in the puo channel, measurements have been performed at ELSA, JLab, 
and GRAAL to obtain high quality data at energies from production thresholds 
to W ~ 2.5 GeV. In this low energy region, the meson-exchange mechanism plays 
an important role and must be treated correctly for extracting N* resonance pa- 
rameters from the data. This is illustrated in Fig. 40. We see that the diffractive 
Pomeron-exchange (dash-dotted curve) becomes negligible at energies near u> photo- 
production threshold. The s- and u-channel nucleon terms, and 7r and r\ exchanges 
can account for the main part of the total cross section. The results shown in Fig. 40 
are from the tree diagram model described in section 4.2 and no N* excitations are 
included. 

y p — > co p 
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Figure 40: 773 — * ujp total cross sections. The theoretical curves are from the model 
of Oh, Titov, and Lee 86 . 




We will here only describe the status of to photoproduction in the resonance 
region( W < 2.5 GeV). The p production will not be discussed since p's width is 
very broad and the coupling of the pN channel to N* states can be meaningfully 
defined only in the analysis involving two pion production channels discussed in 
section 5.4. The <p photoproduction will also not be covered here since <f> meson has 
little, if any, contributions from s-channel resonances, as the ss quark structure of 
the 4> makes N* — > Nip an OZI forbidden decay. 

Quark models that also couple to hadronic channels predict that to photoproduc- 
tion off protons is a promising tool in the search for undiscovered N* states 102 > 101 . 
As in the case of Nrj and K + A, the pu final state, due to the isoscalar nature of 
the u>, is only sensitive to isospin ^ N* resonances. Experimentally, ptu production 
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Figure 41: Differential cross section for "fp — ► pu; in comparison with model predictions of 
Oh, Titov, and Lee 86 . The dotted curves include diffractive production and the uj — > 7r°7 
vertex. The dashed curves are s-channel resonance contributions using the quark model 
predictions of Capstick and Roberts. 



has been measured in both magnetic detectors and in neutral particle detectors. In 
magnetic detectors the pu> channel is usually identified through the u — > 7r + 7r~7r° 
decay. This channel has a 89% branching ratio. Detectors with large acceptance for 
the detection of photons allows to use the to — > tt°-j channel with an 8.5% branching 
ratio. 

The low energy to photoproduction data have been obtained ELSA, JLab, and 
GRAAL. In Fig. 41 we show preliminary differential cross sections from CLAS in 
comparison with predictions using the model of Oh, Titov, and Lee 86 . The model 
contains contributions from diffractive production, ir° exchange, and s-channel N* 
contributions with the photocouplings from the constituent quark model of Cap- 
stick 101 , and the N* — > pui couplings predicted in the model by Capstick and 
Roberts 102 . At high W and forward angles the cross section is completely domi- 
nated by the t-channel processes, i.e. diffractive and pion-exchange contributions. 
Resonance contributions are evident at larger angles, and they seem to play an im- 
portant role in the entire mass range covered by the data. The quark model perhaps 
underestimates the resonance contributions. 

The theoretical models for investigating low energy ui production are still in a 
developing stage. Most of the calculations, such as those displayed in Fig. 40 and 
41, are based on tree-diagrams models, It has been recognized that coupled- channel 
effects must be accounted for before the model can be used reliably to extract 
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Figure 42: Photon Asymmetry of ui photoproduction at 1.125 GeV. The data are 
from GRAAL. The dotted curve is from a tree-diagram calculation. The solid curve 
includes the one-loop coupled-channel effect calculated in Ref.87. 

resonance parameters from the data. The importance of the coupled-channel effects 
on u> photoproduction has been demonstrated in a one-loop calculation 87 based on 
the dynamical coupled-channel formulation Eqs.(52)-(59). As illustrated in Fig. 42, 
the photon asymmetry at E 1 = 1.125 GeV can be changed drastically if the coupling 
with the 7T./V channel is included in the calculation. The K-matrix coupled-channel 
model of the Giessen group 60 , based on Eq.(71), has also been used to investigate the 
data from GRAAL. More efforts are needed to improve these theoretical approaches. 

5.6. Comments on the search for "missing" baryon resonanes 

The search for all baryon states predicted by the SU (6)® 0(3) constituent quark 
model is without question of the highest importance for the field. This model has 
the largest number of excitation degrees of freedom of any quark model based on 
constituent quarks and must form the basis for this search. Hints of possible new 
states, even claims of discovery have been presented in the analyses of single channel 
processes 84 > 161 484, 185,202, 207 ag we u as m a coupled-channel analysis 60 . 

The analysis of the Giessen group is currently the most extensive in searching for 
new states. They employ a coupled-channel K matrix model, as described in section 
4.4.2, and include all available pion and real photon induced reaction channels, 
jN, irN, Nittt, Nr], KA, KY,, Nuo. This analysis finds evidence for several new 
states. While there are clear indications of new resonances near 1900 MeV in some 
data sets, various analyses do not allow to draw definite conclusions on the partial 
waves that are needed to explain the data. The Giessen model may be the most 
promising approach in the search for new states, however, the large amount of data 
fitted simultaneously with many fit parameters involved makes it difficult to assess 
the systematic uncertainties in the fit. The dynamical model approaches should 
be complemented by more experiment-oriented techniques as the one described 
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in section 5.4. Here the binning of data and the evaluation of one-dimensional 
projections of the multi-dimensional parameter space is replaced by an event-by- 
event analysis in specific partial waves that retains the correlations in the data for all 
variables. Theoretical efforts are needed to provide realistic background amplitudes 
that could strengthen the reliability of these techniques. 

6. Concluding remarks and outlook 

In the past few years, we have witnessed very significant progress in the study 
of N* physics. We now have fairly extensive data for ir, 77, K, ui and irir produc- 
tion channels. Much more data will soon be available. The theoretical models for 
interpreting these new data and/or extracting the N* parameters have also been 
developed accordingly. 

^From the analyses of the single pion data in the A region, quantitative in- 
formation about the -fN — > A transition form factors have been obtained. With 
the development of dynamical reaction models, the role of pion cloud effects in de- 
termining the A excitation has been identified as the source of the long-standing 
discrepancy between the data and the constituent quark model predictions. More- 
over, the <5 2 -dependence of the — > A from factors has also been determined up 
to about Q 2 ~ 6 (GeV/c) 2 and found to be in good agreement with the predic- 
tions from a dynamical model. The extracted Ml, E2 and C2 7iV — > A transition 
form factors should be considered along with the proton and neutron form fac- 
tors as benchmark data for testing various hadron models as well as Lattice QCD 
calculations. 

The combined analyses of the ir and 77 production data had led to a rather quan- 
titative, perhaps nearly model independent, determination of several N* parameters 
in the second resonance region. However, a correct interpretation of the extracted 
TV* parameters in terms of the current hadron model predictions requires a rigorous 
investigation of the dynamical coupled-channel effects which are not included in 
the employed amplitude analyses based on either the K-matrix isobar model or the 
dispersion relations approach. 

The analyses of the K, ittt and w channels are still in the developing stage. 
So far, most of the analyses are based on the tree-diagram models with the isobar 
parameterization for the N* excitations. The final state interactions, as required by 
the unitarity condition, are either neglected completely or calculated perturbatively 
using effective Lagrangians. The coupled-channel K-matrix effective Lagrangian 
model, pioneered by the Giessen Group, looks very promising for extracting the N* 
parameters from a combined analysis of all channels. But much works are needed 
to reduce the uncertainties in their non-resonant parameters and to account for the 
ttttN unitarity condition. For a rigorous interpretation of the extracted resonance 
parameters in terms of predictions from hadron models or Lattice QCD calculations, 
the analyses based on the dynamical coupled-channel model, as given in section 4, 
are indispensable. Progress is being made in this direction. 

In the search for new baryon states, progress is being made in developing partial 
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wave analysis procedures that make full use of correlations in the multi-dimensional 
phase space presented in the complex final states. This effort must be supported 
by the development of the theory to obtain improved descriptions of background 
contributions to specific partial waves. 

To end, we mention that the progress we have made in the past decade resulted 
from rather close collaborations between experimentalists and theorists. With much 
more complex data to be analyzed and interpreted, such collaborations must be 
continued and extended in order to bring the study of N* physics to a complete 
success. 
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